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ABSTRACT OF THE DISSERTATION 



Effects of Quenched Randomness on Classical and 
Quantum Phase Transitions 

by Rafael L. Greenblatt 
Dissertation Director: Joel L. Lebowitz 

This dissertation describes the effect of quenched randomness on first order phase tran- 
sitions in lattice systems, classical and quantum. It is proven that a large class of 
quantum lattice systems in low dimension (d < 2 or, with suitable continuous sym- 
metry, d < A) cannot exhibit first-order phase transitions in the presence of suitable 
("direct") quenched disorder. 



ii 



Acknowledgements 



Nothing is ever truly the work of a single person. Much of the content of this disserta- 
tion was shaped in the course of conversations with colleagues and many distinguished 
scientists, among them Royce Zia, Moshe Schechter, Enza Orlandi, Vieri Mastropi- 
etro, Tobias Kuna, Michael Kiessling, Sheldon Goldstein, Alessandro Giuliani, Gio- 
vanni Gallavotti, Diana David-Rus, Eric Carlen, Gabe Bouch, and especially Michael 
Aizenman. 

Anything else achieved here would have been impossible without the support and 
encouragement of my friends; besides those I have already listed above, I owe special 
thanks to Joe Cleffie, Anand Gopal, Alicia Graham, Sarah Grey, and Pankaj Mehta. I 
am also grateful to my brothers and sisters in the Rutgers Council of AAUP Chapters- 
AFT for all they have done to improve the working lives of graduate employees like 
myself. 

Ron Ransome, the director of the physics graduate program for most of my time 
here, has been extremely helpful at many points, particularly in making it possible for 
me to resume my studies after an interruption. 

Above all, I thank my advisor and mentor, Joel Lebowitz, who has helped me (like so 
many other students) in more ways than I can express here, not least by the remarkable 
example he has set for all of us. 



iii 



Dedication 



To my first math teacher and my first physics teacher: 
my mother, Susan Selvig 
and my father, Richard Greenhlatt 



iv 



Table of Contents 



Abstract 

Acknowledgements 

Dedication 

List of Figures 

1. Introduction 

1.1. The rounding effect for classical systems 

1.1.1. Ising models 

Nearest neighbor Ising chain at zero temperature 

Long range interactions in one dimension 

Higher dimensions 

1.1.2. The 3 dimensional XY model 

1.2. The rounding effect for quantum systems 

1.2.1. Transverse field Ising models: direct and orthogonal randomness 

1.2.2. The quantum Ashkin- Teller chain: an exception? 

2. Proof of the rounding effect: overview and preliminaries 

2.1. Notation and systems under consideration 

2.1.1. Systems with continuous symmetries 

2.2. Thermodynamic limit and notions of long range order 

2.3. Statement of main results 



3. A nonlinear central limit theorem 

3.1. Background 

3.2. Definitions 



3.3. The proposition 

3.4. Proof of Proposition 13. 3lT] 

3.5. Bounds on 

4. Free energy fluctuations 

4.1. Definition of Gl 

4.2. Proof of proposition 14. 1 . 1] - finite temperature 

4.3. Proof of proposition 14. 1 . 1] for absolutely continuous distributions of ?7 . . 

4.4. Proof of Proposition 12. 3T] 

4.5. Systems with continuous symmetry: Proof of Proposition 12.3^ 

5. Conclusion 

Appendix A. Methods for numerical studies of random field spin systems 
A.l. The maximum flow representation of the Ising model ground state . . . 
A. 2. Monte Carlo methods for XY and clock models 

A. 2.1. The lookup table algorithm for the Clock model 

A. 2. 2. A modified Ziggurat algorithm for the XY model 

Appendix B. Rounding of First Order Transitions in Low-Dimensional 
Quantum Systems with Quenched Disorder 

Appendix C. Some concepts and results in mathematical probability . 

C.l. a- algebras and measures 

C.2. LP norms and spaces; convergence of measurable functions 

C. 3. Random variables, expectations, conditional expectations 

Appendix D. Product Measure Steady States of Generalized Zero Range 
Processes 

D. l. Introduction 

D.2. Factorizability in the Mass Transport Process 

D.3. Reverse processes 



3C 



33 



37 



39 



3S 



41 



44 



46 



48 



53 



54 



54 



56 



56 



58 



61 



3 



vi 



D.4. Factorizability in Generalized Zero Range Processes 

D.5. GZRPs on infinite lattices 

D.6. Conclusion 

References 

Vita 



vii 



List of Figures 



1.1. Plot of ground state magnetization as a function of mean magnetic field 
/iQ for a random field Ising chain with magnetic field distribution ^6f^_jj + 
i^+j^, J = 1,H = 0.42 y 

1.2. A two-dimensional example of a dichotomous random field configuration 
leading to "Perestroika" in the ground state. The bold lines divide sepa- 
rate regions where the spins are always -|-, indeterminate, and always — 

in any ground states 

A.l. Example of the maximum flow graph used to derive Equation (|1.1.5p . Q 
A. 2. The ziggurat algorithm: a probability density, bounding boxes of equal 

volume l59 

A. 3. The modified ziggurat algorithm: two probability densities, bounding 

boxes of equal volume. The probability densities shown are those of 

Equation (lA.2.3p . with k = 1 and k = 1.5, rescaled for pfc(O) = 1 l60 



viii 



1 



Chapter 1 
Introduction 



One of the basic techniques of condensed matter physics is the effective description 
of sohds as a combination of a static portion and a rapidly- moving part: in a simple 
description of metals, the nuclei and tightly-bound electrons remain fixed and form an 
effective potential background for the conduction electrons. Although the simplest de- 
scription involves a uniform or periodic background (i.e. a perfect crystal), this is hardly 
a natural assumption: completely pure samples are anything but common or easily pre- 
pared. There are many situations in which disorder has only a minor effect, but there 
are cases of fundamental importance where this is not the case. The best-established 
illustration is the description of electrical conduction in metals: with the application 
of quantum mechanics in this context it became clear that irregularly-placed scatter- 
ers were necessary to account for finite conductivity. Anderson localization provides 
a further way in which disorder produces a qualitative difference in the behavior of a 
physical system. 

These examples concern transport properties, which are harder to fit into a compre- 
hensive framework than equilibrium properties. The core of the work described here is 
a similarly qualitative effect at the level of equilibrium thermodynamics, the rounding 
effect predicted by Imry and Ma [Ij and described in detail in the next section. 

It is misleading in a way to talk about equilibrium in this context. As noted already 
in the paper which introduced the mathematical framework now known as quenched 
randomness [2] , it is important to consider a background which is a metastable configu- 
ration, and is not typical of the equilibrium state of the full system. Although disorder 
is still present in systems which are genuinely in equilibrium (this is what is known as 
annealed disorder), annealed systems cannot exhibit behavior which is fundamentally 
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different from that of ordered systems. 

There is still much that remains to be understood about classical models of quenched 
randomness, but the situation for quantum systems has been even more obscure. The 
main results presented in this dissertation, Propositions 12 . 3 .T] and ^.3.2 \ are an extension 
of the proof of the Imry-Ma rounding effect to quantum systems. The present chapter 
will discuss the previous state of understanding and attempt to provide context for the 
result. Chapters [2] to H] comprise the proof of these results. Chapter [2] describes the 
formalism used, establishes several preliminary results, and states the main proposi- 
tions to be established. Chapter [3] contains a proof of a nonlinear central limit theorem 
(based on an earlier result of Aizenman and Wehr [3]) which may be of some inde- 
pendent interest. Chapter H] completes the proof of the main results with an analysis 
of the free energy effects of the quenched randomness. This work was announced in 
a publication by the author with M. Aizenman and J. L. Lebowitz, which provides a 
summary of the argument as is therefore attached as [4J. Additionally, Appendix [Q 
reviews some probabilistic terminology and results used in the previous chapters which 
may be unfamiliar to some readers. 

Appendix [D] (written with J.L. Lebowitz, and published as [5]) describes earlier 
work by the author, with J. L. Lebowitz, on nonequilibrium stochastic dynamics. 

1.1 The rounding effect for classical systems 

A 1975 paper by Imry and Ma contains an important insight into phase transitions 
in disordered systems based on an analysis of the energy of the ordered phase, using 
considerations similar to those applied more rigorously in Peierls' proof of long range 
order in the Ising model [6] and later in Pfister's proof of the Mermin- Wagner theorem 
for classical systems [7j. The context is 0{N) models, that is lattice models where 
configurations consist of a specification of an A^-dimensional unit vector (a classical 
spin) (7x at each site x, with equilibrium states determined by the Hamiltonian 



Since the paramagnetic phase of this system has higher entropy, for the ordered 
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phase to be stable requires that the energy cost involved in forming a domain where 
the spins inside are aligned in a different direction from those outside grow with the 
size of such a domain. In the Ising {N = 1) case the cost for a domain of diameter L is 
of order L'^~^, and in continuous (A^ > 2) versions spin- wave analysis [8j suggests a cost 
on the order of L'^^^. In the absence of a random field this suggests (correctly) that 
ferromagnetism does not exist in these systems at finite temperature for d = 1 and d <2 
respectively. Since ferromagnetism appears at any higher dimensioj^, we may speculate 
that it is sufficient for the energy cost to grow with L, a contention that is supported by 
estimates of the number of genuinely independent contours of given size [9l|T0]. If the 
random field has typical strength H and we neglect correlations between the field at 
different locations, then the total random field in a domain of volume L'^ will typically 
have magnitude HL'^^'^. Then when d < 2 for Ising models and d < 4 for continuous 
models, and given any direction, there will be a large number of large domains for 
which flipping into that direction is energetically favored. On this basis, Imry and Ma 
predicted that there would be no long range order at low temperature for the random 
field Ising model in two dimensions and for similar continuous models in d < 4. They 
also suggested that ferromagnetism would persist in higher dimensions. 

Another way of looking at ferromagnetic order in this system first order 

transition, where the equilibrium magnetization (a) changes discontinuously as the 
external field h is changed through zero. The disappearance of ferromagnetic order 
corresponds to a "rounding" of this discontinuity, leaving a continuous transition. We 
shall see that this "rounding effect" occurs in a large number of systems in the presence 
of quenched randomness. 

In 1976 Aharony, Imry and Ma established a detailed connection between random 
field 0{N) models with continuous spin in 4 < d < 6 dimensions and the field- free 
versions in d — 2, finding an exact correspondence between the most divergent Feynman 
diagrams of all orders for the two models [llj; among other things this provided strong 
support (which had previously been lacking) for the prediction that the random field 



^With the possible exception of the symmetric quantum case, where a ferromagnetic phase has yet 
to be rigorously shown to exist in three dimensions. 
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models had ferromagnetic order for d > 4. However by expressing the Lagrangian of 
the model in a supersymmetric form, Parisi and Som'las were able to extend this per- 
turbative correspondence to all dimensions and to n = 1, suggesting that the random 
field Ising model was not, in fact, ferromagnetic in three dimensions but only for four 
dimensions or more |12j . or in other words that its lower critical dimension di was 3. A 
number of attempts to study the formation of domain walls more carefully than Imry 
and Ma seemed at first to agree on = 3 [131 ttH but before long other domain- 
wall studies appeared to return to = 2 [UllTT], along with other theoretical |18j and 
experimental work [19j; in particular Chalker [9] and Fisher, Frohlich and Spencer [10] 
provided strong (but not conclusive) arguments for di = 2 based on a rigorous treat- 
ment of the "no contours within contours" approximation. However further arguments 
emerged for = 3 |20j . and the debate was only resolved with rigorous proofs of long 
range order for the 3 dimensional random field Ising model by Imbrie [211 [22] (for 
zero temperature) and Bricmont and Kupiainen \12>\ 124] (for low temperature), based 
on intricate examinations of the scaling behavior of the contour representations of the 
model. 

This did not yet completely vindicate Imry and Ma's argument; this was done 
by Aizenman and Wehr, who proved that first order transitions could not exist for a 
large variety of classical systems in the presence of disorder [251 IS]- They were able 
to do this by first constructing a suitable description of the equilibrium states of the 
infinite system (metastates) , which allowed the construction of a quantity describing 
the free energy fluctuations due to the random term in the Hamiltonian. The estimates 
of domain energies in the Imry-Ma argument correspond to rigorous bounds on this 
quantity, and by examining only hypercubic domains it is possible to show that a first 
order transition would cause a contradiction between these bounds in the dimensions 
which Imry-Ma predicted a rounding effect, that is always in d < 2, and for systems 
with continuous symmetries d < 4. 

The precise conditions are somewhat cumbersome to state precisely. They are ex- 
actly the same as those of Propositions 12.3.11 and 12.3.21 below, so for the moment I will 
confine myself to some general remarks. The main one is on the relationship between 
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the quenched disorder and the order parameter. If the phase transition or long range 
order under examination is described by averages of some local quantity Kx, then a 
rounding effect can be expected only when the Hamiltonian can be written in the form 



following Hui and Berker |26] , we can refer to this as "direct randomness" . In the Ising 
model, a random field is direct with respect to the spins, and so an arbitrarily weak 
random field eliminates ferromagnetism in two dimensions; bond randomness is not 
direct, so it does not (at least not when it is sufficiently weak). It should be noted that 
this notion is relative to a particular phase transition: by way of illustration we may 
consider the random-bond Potts model in two dimensions. In the nonrandom version 
of this model with sufficiently many colors, the order-disorder transition involves a 
nonzero latent heat, which in this case means that the equilibrium bond energy density 
is discontinuous with respect to the bond strength. This is a first order transition for 
which bond randomness is direct, and therefore the latent heat vanishes whenever it is 
present. On the other hand this randomness does not couple to the color, and so as in 
the Ising model long range order remains. 

1.1.1 Ising models 

More can be said about the random field Ising model (henceforth RFIM) by bringing 
a variety of specialized techniques to bear, leading to more insight into the scope and 
significance of the rounding effect. In particular, some insight can perhaps be obtained 
into what replaces the ferromagnetic phase in this situation, and into the more subtle 
effects of the random field in 3 dimensions. 

Besides computational efficiency, one major issue in simulating disordered systems 
is the presence of additional finite size effects, especially at the lower critical dimension. 
To make this clearer, let us revisit the Imry-Ma [Ij analysis of the random field Ising 
model with random fields of typical strength H. Flipping a typical domain of linear 
size L will involve a bond energy of the order JL'^^^ and field energy on the order of 
in one dimension, the field energy will dominate once length scales on the order 




(1.1.2) 



X 
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of (J/H)'^ come into play, but when smaller systems are analyzed they will appear to be 
ferromagnetic. In two dimensions the competing energies are both proportional to L; 
when H is small compared to J the rounding effect occurs only because of fluctuations 
in the random field, which makes its effect stronger in particular regions. This can be 
studied by means of extreme value statistics, and this approach [17J gives a breakup 
length scale on the order of 



with a constant A of order 1. This has been backed up by numerical studies, which 
found A = 2.1 lb .2 for a Gaussian distribution of the random fields and 1.9 it .2 for a 
bimodal distribution [27]. For weak values of the random field, this distance can easily 
be hundreds or even thousands of sites - nowhere near macroscopic, but potentially 
very difficult to reach in simulations. 

Nearest neighbor Ising chain at zero temperature 

One requirement of the Aizenman-Wehr proof of the rounding effect are assumptions 
which must be made on the distribution of the random parameter in certain contexts. 
Some limitations may be purely technical (see Section 13.51 below), but not all. We 
can see this thanks to studies of the one-dimensional Ising model by Bleher et. al. |28] . 
Examining the case of a "dichotomous" random field, i.e. one taking only the two values 
zizH and those with equal probability, they found that the ground state configuration 
of the spin at any site x could be deduced from the random field in some finite but 
undetermined neighborhood as follows. 

Let us write the Hamiltonian of the system as 



Lt = exp [A{J/Hf] 



(1.1.3) 




(1.1.4) 



X 



X 
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We can recursively define two position-dependent functions of the random fields by 

Ux-l+hx, \Ux-l + hx^i\ < J 

J, Ur,_i + > J (1.1.5) 

-J, Ux-l + hx-1 < - J 

Vx+l + hx, \Vx+l + hx+i\ < J 

J, Ux+i + hx+i>J ■ (1.1 -6) 

— J, Ux+l + hx+l<—J 

Ux (respectively Vx) can be thought of as representing the effect of x's neighbors to the 
left (resp. right) on flipping it out of the ground state - the lowest-energy flip may involve 
a number of sites depending on the magnetic field they experience. These quantities 
always exist, and are almost always uniquely specified since there will eventually be a 
large block of sites where all the magnetic fields point in the same direction. Any site 
X for which Ux + Vx + hx is positive (resp. negative) will necessarily have ax = 1 (resp. 
— 1) in any ground state, and Ux + Vx + hx = there will be ground states with ax = 1 
and ax = —1. All this is proven (for ijx = ^H) in [28], but readers should be able to 
convince themselves by considering the minimum energy cost of flipping a block of sites 
containing x out of the resulting configuration; and Appendix |X] contains a derivation 
for arbitrary fields. 

It is not difficult to numerically estimate the probability distribution of uq from the 
recursion relationship (jl.l.Sp (the distribution of vq is identical and independent), and 
from this calculate the average value of the ground state magnetization. Figure 11.11 
shows a plot resulting from such a calculation 

As is apparent from Figure 11.11 the magnetization in the presence of a dichotomous 
random field has a number of discontinuities, in fact an infinite number occurring 
wherever h/H is rational. It is interesting to note that these first order transitions do not 
correspond to any long range order: there are a finite density of isolated regions which 
can be flipped independently with no change in energy, resulting in a finite residual 
entropy; this situation was called "Perestroika" when first described in 1989 [25] . 

Nonetheless, this illustrates that one of the restrictions on the proof of the rounding 



-1-0 - 



Figure 1.1: Plot of ground state magnetization as a function of mean magnetic field ho 
for a random field Ising chain with magnetic field distribution ^6f^_jj + ^dj^_^_jj, J = 1, 



effect by Aizenman and Wehr [3], the requirement of an absolutely continuous distri- 
bution of the random field at zero temperature, is indeed necessary. 

These first order transitions do not appear for absolutely continuous distributions 
of the random field. The curves do not, however, appear to be always analytic, and the 
character of the singularities (that is, the order of phase transitions present) appears to 
depend on the corresponding properties of the random field distribution in a way that 
remains to be investigated more carefully. 

Long range interactions in one dimension 

The one dimensional Ising model can exhibit long range order at finite temperature if 
interactions are sufficiently long range [30j. Let us consider the variant of the RFIM 
with the following Hamiltonian: 



H = 0.42 . 




(1.1.7) 
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Figure 1.2: A two-dimensional example of a dichotomous random field configuration 
leading to "Perestroika" in the ground state. The bold lines divide separate regions 
where the spins are always +, indeterminate, and always — in any ground states. 

The bond energy associated with flipping a block of L spins is on the order of 
so the Imry-Ma argument indicates that rounding should occur for a > 3/2. This was 
confirmed by Aizenman and Wehr [3] , but the question of what happens for even longer 
ranged interactions remained unanswered until recent work by Cassandro, Orlandi and 
Picco [31], who showed that long range order persists in the presence of weak random 
fields for 1 — ln(3/2) < a < 3/2. This means that the estimate provided by the Imry-Ma 
argument is also sharp in this respect, and suggests that the restrictions on long range 
interactions used below may in some sense be sharp as well. 

Higher dimensions 

There is a complication in the ground state behavior of the RFIM which has not been 
well studied: "Perestroika" (see p. [7]) occurs in the ground state for dichotomous ran- 
dom fields in all finite dimensions, because there will be a finite density of regions where 
the magnetic field has a pattern like that shown in Figure [L2l In two dimensions, this 
means that a dichotomous RFIM will exhibit a first order transition at zero tempera- 
ture, which has probability zero [5] either for the same model at finite temperature or for 
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any absolutely continuous distribution of the random fields. In three dimensions there 
is also a more conventional degeneracy in the ground state due to long range order [22] , 
but numerical studies still indicate that dichotomous and absolutely continuous distri- 
butions show strikingly different behavior, not even lying in the same universality class 
at zero temperature |32[ 133] . If Perestroika is the main cause of the difference between 
the two cases then it should disappear at finite temperature, with all low temperature 
systems behaving like the zero temperature system with an absolutely continuous field 
distribution. This appears to be supported by comparing finite-temperature Monte 
Carlo studies with dichotomous [M] and Gaussian [35j random fields, as well as com- 
paring the latter to ground state studies with Gaussian fields [35]. It is possible that 
renormalization studies which include parameters differentiating between the different 
distributions could shed light on the situation; then the scenario described above would 
involve additional dichotomous-field fixed points, all with this new parameter as an 
unstable direction. 

Whatever the details, it is clear that different types of random field distribution can 
result in profoundly different types of behavior, especially in the ground state. 

1.1.2 The 3 dimensional XY model 

A claim has arisen recently that the prediction of a rounding effect for three and four di- 
mensional systems with continuous symmetry (where dimensional reduction, the Imry- 
Ma argument, and the proof of Aizenman and Wehr are all in agreement) is either 
incorrect or misunderstood. The controversy has been specifically about what is prob- 
ably the simplest such model, the three dimensional random field XY model described 
by the Hamiltonian 

m = -j Yl • - ^ 4 • (1.1.8) 

<x,y> X 

where ax are unit vectors in M?, and hx are i.i.d. random vectors in M^; this is the 
N = 2 case of the 0{N) model discussed above. In the following discussion, we can 
assume that hx are chosen uniformly from some circle of specified radius H, as is done 
in most of the numerical studies we will discuss. 
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In 2007, Fisch |36J published the results of Monte Carlo simulations on the random 
field clock model, which shares the Hamiltonian (jl.l.Sp . but where a and h are now 
restricted to a set of q evenly-spaced directions (in [36] g = 12 is used). Although this 
model does not have the continuous U{1) symmetry of the XY model, there is certainly 
a relationship between the properties of the two models [Ml EZl EH] , and the clock model 
can be simulated very efficiently. 

We can be more concrete in considering the nonrandom [H = 0) model in two dimen- 
sions. Here we know that the XY model has no long range order (i.e. no ferromagnetic 
phase) at finite temperature thanks to the Mermin- Wagner theorem [39 1 140 1 ^T] . The 
clock model, on the other hand, has a finite number of ground states, which are related 
by a symmetry group with minimum interface energy of (1 — cos 27r/g) J per bond, and 
so by Pirogov-Sinai theory [411 H2] have ferromagnetic long range order for sufficiently 
small temperatures. 

An early Monte Carlo study of two dimensional clock models [38] noted that outside 
the ferromagnetic phase the clock model behaved similarly to the XY model, in particu- 
lar showing evidence of a Kosterlitz-Thouless phase of quasi- long-range order. It is not 
surprising that the relationship between the two systems should depend significantly 
on the temperature, since the relatively high energy excitations will be similar between 
the two systems. 

It was claimed in [36] that the temperatures under consideration were high enough 
that differences between the clock and XY models would not come into play, but there 
are reasons to doubt that this is a reasonable line of argument. Ferromagnetic order 
can always be disrupted by a system's lowest energy excitations; their effects only 
become less relevant with increasing temperature insofar as they are overwhelmed by 
other, more entropically favorable, excitations. If the QLRO phase begins at a nonzero 
temperature, it is because it is only then that the associated modes begin to play a 
dominant role, and it is at exactly this temperature that the absence of sufficiently 
low-energy excitations in the clock model makes itself felt. 

It is also worth noting that some studies of the clock model |43] involve some sites 
with zero magnetic field; although ferromagnetism is still ruled out at finite temperature 
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(see Section [2.1.11 below) . this does open the possibiHty of a ferromagnetic ground state 
or some other phenomenon similar to perestroika which could have subtle effects at 
finite temperature. 

It appears to be possible to virtually eliminate discretization effects with another 
scheme which at least provides substantial improvements over rejection sampling. The 
idea, described in Appendix lAl is based on the Ziggurat algorithm [S], a method which 
has proven to be highly efficient in sampling the normal distribution |45j . 

Preliminary tests of this method have been very promising. Using a C++ program 
on a desktop computer with a 3.2 GHz Pentium 4 processor and 2 GB RAM, I have 
been able to achieve an update rate of 8.6 x 10^ to 2.1 x 10^ sites per seconcj^ on a 
parameter range k E [0, 14], compared to 6.8 x 10^ to 1.3 x 10^ for a comparable lookup- 
table implementation of the 12-state clock model. It seems very likely, then, that it will 
be possible in the near future to conduct a detailed study comparing the two models. 

1.2 The rounding effect for quantum systems 

1.2.1 Transverse field Ising models: direct and orthogonal random- 
ness 

The simplest quantum lattice spin systeno is the transverse field Ising model, defined 
by the Hamiltonian 

^ = - ^ JxyO'3,x(^3,y " ^ ^xC^l^x - ^ /ia;0"3,x, (1-2-1) 

where ai^x denotes the i component of a ^-spin at site x. The properties of this system 
(and a number of variants) are relatively well known, in large part due to the fact 
that its path integral representation is the continuum limit of an Ising model with 
an additional dimension (sometimes called the "space-time Ising model") |49[ I50j. Its 



The efficiency depends on tiie system size, which may indicate that further optimization is possible. 

^That is, the simplest lattice spin system which involves nontrivial commutation relationships. The 
Ising model, for example, has no classical dynamics, and is in a certain trivial sense a quantum system 
- the DLR conditions [46] which define its equilibrium state are equivalent to the quantum KMS 
conditions defined by Heisenberg evolution [47). This is unlike ofT-lattice systems where the classical 
and quantum KMS conditions do not coincide [48] . 
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behavior is consequently very close to that of a classical Ising model in many respects; 
the nonrandom version is ferromagnetic when A is small, for example. One difference is 
that the phase diagram of the system at zero temperature is more complicated than that 
of its classical counterpart; the system has a ferromagnetic-paramagnetic transition at 
zero temperature at a critical value of A, which provides a paradigmatic example of a 
quantum critical point |51j . 

Among the ways of introducing quenched randomness to this system, the most 
straightforward are to add randomness in the transverse field A^ or the longitudinal 
field hx- We can hardly expect quantum effects to be very striking in this system, but 
it is still worth clarifying where it fits into the picture I have been discussing. 

The random transverse field case is a good example of orthogonal randomness, in 
that ferromagnetic order remains as long as the transverse field is not too strong [52j ; the 
nature of the ferromagnetic-paramagnetic transition can be changed significantly [53], 
but the details of this are beyond the scope of the present work. 

A random longitudinal field, on the other hand, couples to the magnetization, and 
so should be direct randomness. It has been expected [Mj that the outcome should be 
similar to the (classical) random field Ising model, to which it reduces for A^ = 0. 

1.2.2 The quantum Ashkin- Teller chain: an exception? 

The possibility of additional complications in the quantum case have been raised in the 
context of the quantum A^-color Ashkin- Teller model. In this system, each lattice site 
X contains ^-spins, described by operators a'f^^ for the ith component of the a spin 
at site X. In one dimension, the Hamiltonian is given by 

N 

' (1-2.2) 

a<f} X 

To begin with, we examine the nonrandom version of the system, where Jx = J and hx = 
h are constant. At sufficiently low temperature (zero temperature in one dimension) 
and when J is large compared to /i, the system is in an ordered "Baxter phase", with 
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the spins of each color exhibiting long range order, with no simple correlation between 
the different colors, while for large h the system is paramagnetic. For > 3 and e > 0, 
the transition between these states is of first order [55], characterized for example by a 
discontinuity in (^cr^^'^ (which is independent of x and a). Randomness in hx is clearly 
direct with respect to this transition, and so is randomness in Jx - the two can be shown 
to be equivalent by a duality transformation [55]. Therefore we should expect a system 
with such randomness to round the first order transition (at least provided it has an 
absolutely continuous distribution, cf. Section [1.1.11 above). 

A renormalization group analysis by Goswami, Schwab and Chakravarty [55] sug- 
gested that this might not be the case. They found that when e was below a certain 
nonzero value ec{N), the flow of the system was similar to that of the random transverse- 
field Ising model and that there was no first order transition; however above this value 
their scaling analysis broke down in a way that led them to suggest that a first order 
transition might persist. As we shall see this can be rigorously ruled out, but we are 
not yet in a position to say exactly what is happening. 
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Chapter 2 

Proof of the rounding effect: overview and preUminaries 

We now embark on the proof of the rounding effect. The basic framework of the 
argument is the same as [3J, which in turn uses reasoning based on that of Imry and 
Ma [Ij. One constructs a random variable Gl which represents the free energy effect 
of the random field on a scale L. We then show that it has a strict upper bound of the 
form 

\Gl\<CL'^'^ + C'L'^/'^ (2.0.1) 

or in more restricted cases 

\Gl\< CL'^-^ + C'L'^/^. (2.0.2) 

At the same time, we show that when the system is at a first order transition, it has 
asymptotic fluctuations described by a normal distribution, 

Gl f« AA(0,L'^/2) (2.0.3) 

on the scale L*^/^, which means that it will violate the above bounds in sufficiently low 
dimension. 

The behavior indicated in Equation (|2.U.3|) is akin to a central limit theorem, but 
instead of a sum of random variables it concerns a suitably continuous function of a 
large number of random variables. In Chapter [3] we present a suitable nonlinear central 
limit theorem. This result is a slight modification of one found in [3j. Although the 
result is phrased in what we hope will be a more useful form for some readers, the proof 
is substantially the same, apart from a correction due to Bovier |56j . 

The upper bound (j2.0.ip is quite easy to show for finite systems, however it is 
not trivial to show that an infinite-system limit exists. This problem was resolved for 
classical systems by defining as expectation values with respect to metastates, which 
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are random probability measures related to the random Gibbs states of a disordered 
classical system [56]. The notion of metastate has been generalized to one suitable to 
quantum systems (that is, one based on the operator analysis notion of KMS state rather 
than the measure-theoretical notion of Gibbs state) by Barreto and Fidaleo \57\ I58j . 
but while this is promising for many other problems in disordered systems it is of little 
use to us. Instead, we have formulated an argument which remains almost exclusively 
at a thermodynamic level. This has the additional merit of producing a proof which is 
considerably more accessible from both a physical and a mathematical point of view. 
Chapter m begins with the construction of an object satisfying the upper bound (j2.0.ip 
and the conditions of the nonlinear central limit theorem proven in Chapter [3l completes 
the proof of our first main result, and then provides the additional estimates needed 
to obtain the bound (j2.0.2p under suitable conditions and obtain a stronger result for 
systems with continuous symmetry. 

Before embarking on the proof, we establish definitions and a number of preliminary 
results which establish the context, and state our two main results. 

2.1 Notation and systems under consideration 

We consider systems on a lattice (we take this to be the simple cubic lattice Z"^ for sim- 
plicity, but many other cases can be reduced to this), where the possible configurations 
of each site are described by a finite-dimensional Hilbert space, with time evolution 
affected by a static background described by means of its statistical properties. 

To make this more mathematically precise, we suppose that we are given a dimen- 

^0- We introduce a copy Ax of this 
algebra for each lattice site x S Z"^, and take everything which can be obtained by 
tensor products, sums, and limits: this is the quasi-local C*-algebra A [59l|60lll7], and 
we will take the conventional point of view that this allows us to describe all physical 
observables. We let F be the finite subsets of Z*^, and for any A € F we let be the 

C*-algebra is a collection of operators with addition, multiplication, conjugation, which is closed 
under all of these operations (e.g. the product of two operators is another operator in the same algebra) 
and with a norm which defines limits, convergent series, etc. This is a common way of formalizing the 
notion of the set of operators describing a quantum system. [461 1591 l6Uj 
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local C*-algebra on A. 

To specify the background referred to above, we will make use of the following 
concepts: 

Definition 2.1.1. A field (on U^) is a map from Z'^ to the real numbers. The set of 
all fields is denoted by E. 

Definition 2.1.2. A random field (on U^) is a collection of random variables indexed 
by the elements of'L'^. A random field is i. i. d. if the random variables it consists of are 
independently and identically distributed. 

If we consider the space £ to have the cylinder-Borel sigma algebra (the conventional 
choice), then an i.i.d. random field is also a iS-valued random variable (by Kolmogorov's 
extension theorem). 

For Chapter [3l as in many other generalizations of the central limit theorem, we 
need a restriction on the moments of the individual random variables: 

Definition 2.1.3. A random variable X is Lyapunov if there is a 5 > 2 such that 
Av \ X\^ is finite. A random field rj is Lyapunov iff each rj^ is Lyapunov. 

Note that an independent, Lyapunov random field defines an array (by restrictions 
to subsets of U^) which satisfies the usual Lyapunov condition, hence my appropriation 
of that name. 

In what follows, strictly separate symbols will be used to denote random fields and 
specific values, r/ will be a random field (consisting of the individual real random 
variables rjx), while C is a specified (nonrandom) element of 8. It is very convenient to 
have a compact way of referring to the random field within a specified subset of Z"^; to 
do so we use the symbol r]\ to refer to the collection of r/^. with x G A; the meaning 
of expressions like = should be clear. This allows a convention we will use for 
conditional expectations: by expressions of the form 

Av[/(r?)|r/A = CA] (2.1.1) 

we mean a conditional expectation of the random variable /(?]) on the sigma-algebra 
generated by specifications of r/A, understood as a function of C,. 
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Dynamics (and equilibrium states) on such a structure are defined by way of the 
concept of an interaction, basically a rule for assigning Hamiltonians to families of 
systems defined on different finite regions. Formally, a nonrandom interaction is a 
function ^'o : F ^ ^ satisfying *o(-'^) € Ax- 

We wish to consider interactions depending on one or more random fields. For 
the matter at hand, we do not need to talk about arbitrary random interactions; it is 
enough to talk about systems where the Hamiltonian on a finite region F G F with free 
boundary conditions is 

XCT x:Ta,AoCT a=l x-.Ta^AcCT 

where T-c denotes translation by x, and *o is assumed to be translation invariant 
{"^o(TxX) = Tx^o{X)). We define other boundary conditions as follows: 

Definition 2.1.4. A boundary condition is a linear map B : ¥ x A A, {T,A) t-^ 

Bi^(A) satisfying 

1. \\By{A)\\ < \\A\\ for allAeA 

2. Br (A) e At for all A e A 

3. Bt{A) = A for all A ^ At 

4. Br{A) = for all A e Arc 

This is a fairly generous notion of boundary conditions, and in particular includes 
fixed and periodic boundary conditions. We denote the Hamiltonian with boundary 
condition B by 

HrlB- = Yl MMX)) + Y^^ + Cx)Brii^x) + Y Yl ^axBrhax). (2.1.3) 

X xedoT a=lxeda^ 

where daT denotes the set of a; G Z*^ for which TxAa contains members of both F and 
F^. 

This allows us to define partition functions by 

4',b(C,cu) := TreM-PH^lB-), (2.1.4) 
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free energy by 



and Gibbs states by 



FlsiCui) := -^logZ^^siCui) (2.1.5) 

h Tr •e^'^^r.'s" 

To avoid a profusion of subscripts, we omit a label for boundary conditions when 
periodic boundary conditions should be understood; and when an integer L appears 
instead of the finite set T it should be understood to represent the (hyper)cubic subset 
of Tj'^ of side length L approximately centered at the origin, i.e. 



-L + 1/2 L + 1/2""^ 



nZ'^. (2.1.7) 



2 ' 2 

It is helpful to observe that the /3 — )• oo limit of the free energy and Gibbs states (for 
the time being, we consider these limits with all other parameters fixed) exist. Indeed, 
when the ground state is nondegenerate, the free energy converges to the ground state 
energy and the Gibbs state converges to the (unique) ground state. Even in the presence 
of degeneracy, these limits provide an equally useful description of the system, and we 
can establish many results simultaneously for finite and zero temperature by taking 
advantage of this. We will therefore take the free energy and Gibbs states to be defined 
for all /3 € [0, cxd], with the values at /3 = oo being the above limits. 

The free energy, as defined in Equation (j2.1.5p . has the following well-known prop- 
erty we will use repeatedly in what follows: 

Lemma 2.1.5 (|46j). For any Hermitian matrices A,B of the same size, 

jlogTre^ -logTi-e^l < (2.1.8) 

The terms of the interaction connecting a finite region to the rest of the system play 
an important role in the arguments of the present work. We denote these by 

Vl--= Y1 mx)+ J2 (^ + Cx)^x + J^ Yl (2.1.9) 

Note that ||i?r(^i'~)|| ^ ^or all boundary conditions, so bounds on the norm of 

the infinite-system operator above give considerable information about finite systems 
as well. 
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Our main result will be restricted to systems which are short range in the following 



sense: 



Assumption 2.1.6. There are constants < C, C" < oo such that 



Av 



<C{l + \h\)L'^-^ + C'L'^/^ 



(2.1.10) 



This may not be very transparent, so we note the following results which provide 
sufficient conditions under which Assumption 12 . 1 .S] is satisfied. 

Lemma 2.1.7. If t] and v are i.i.d. and mutually independent with finite, then 
there is a constant < ci < cxd such that 



Av 



x-.T^AonVLii^^T^Ao} "=1 x:r^AQnr£^{0,T^Ac} 



(2.1.11) 



Proof. The quantity whose norm is being bounded consists oi N + 1 sums, each with 
no more than 2d\Aa\L'^~'^ terms, each bounded in norm by 1 or □ 



When this holds, it means that Assumption 12 . 1 .6] is satisfied iff the following condi- 
tion on is satisfied: 



w 



0,0 1 



X:XnrL^{0,x} 



(2.1.12) 



This is clearly the case when ^'o is of finite range, but also allows some scope for infinite 
range interactions. A convenient condition [3] is 



Lemma 2.1.8. // 



\dX 



V diamAip-^||*o(^)|| < c'L(2-'^)/2 



X90 
diam X<L 



for all L, then Inequality \ 2. 1.1 2\ is true. 



(2.1.13) 



Proof. By the triangle inequality 



i^rii < 



E ll^o(^)ll; 

X:XnrL^{0,x} 



(2.1.14) 
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and in this sum the terms with diameter L or less contribute, at most, 

^ 2ciL^_idiamX ^2dc'L^/2, (2.1.15) 



diam X<L 



and the remaining portion is bounded by 



E L^^j§l\mX)\\ < L^-' E diamxi^||^'o(X)||<c'L'^/^ (2.1.16) 

diamX<L diamX<L 

Putting the two parts back together we have Inequahtv 1 2 . 1 . 1 2 1 with C = {2d+ l)c'. □ 

For pair interactions, the bound in Lemma [2.1.8l is satisfied in d = 1 for interactions 
decaying like (distance) or faster; a result of Cassandro, Orlandi and Picco [31J 
shows that Proposition 12.3.1] is false for a system with slightly longer range interactions, 
which suggests that Assumption 12.1.6] may in some sense be optimal. This may be of 
some practical interest, since for pair interactions in d = 2 we need the interactions 
to decay strictly faster than (distance)"^ for Lemma [2.1.81 to apply, and inverse cube 
interactions seem to be quite common [6IJ. 

Finally, we give a similar statement which provides some control over the case of 
infinite-range random interactions: 

Lemma 2.1.9. Let v be i.i.d., with Nq, = 00 and 

diamA„|a^„| Av|t;„,o| < cL(2-'^)/2. (2.1.17) 

«>i 

diam Aa <L 



Then 



Av 



^ ^ ^ ^ "^ax lax 

a=l x:T^AanrLt{1>,T^Ac} 



< (2.1.18) 



Proof. The contribution of terms with diamyl^ < L is bounded by 
2d(L + diam ^a)'^^-^ diam Av |i;q,o| 



diam Aa<L 



< d2'^L^-^ diam ^„ | dA^ \Av\vao\< d2'^cL'^'^ , 



(2.1.19) 



diam j4tv<I/ 
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while the remaining terms are bounded by 

L'^\dAa\Av\vao\<L'^'^ diamAa\dAa\Av\vao\<cL'^^^, (2.1.20) 

diam Aci>L diamAa>L 

and the conclusion follows with c' = {d2'^ + l)c. □ 
2.1.1 Systems with continuous symmetries 

Imry and Ma's initial work [1] mainly concerned systems with continuous symmetries. 
In this context the Mermin- Wagner theorem |39i HQ] already precludes long range order 
without randomness in two dimensions, so the rounding effect would be of little conse- 
quence except that it extends to four dimensions, but only so long as the randomness 
preserves the symmetry "on average" in a sense the following passage should make 
clear. 

First, we assume that the single-site algebra contains a subalgebra isomorphic 
to the rotations SO{N) for some N > 2. For each rotation R £ SO{N), let be the 
corresponding element of Ax- We will say that an interaction ^ is invariant iff 

*(x) = ( n ^^'] *(^) ( n ) (2.1.21) 

\xex / \xex / 

for all X G F and aU R G SO{N). 

Intuitively, for a random system to be (stochastically) invariant under rotations, 
the field and the quantity it couples to should both transform as dual representations 
of SO{N). The vector representation is the only case we are aware of which includes 
any cases of intrinsic interest (this case, in particular, includes Heisenberg models in a 
random magnetic field), so we will focus on this. The fields are then elements of £^ , 
or equivalently maps : Z'^ — )• M^, and we let 

Definition 2.1.10. A random vector field is a collection of -valued random variables 
indexed by the elements ofZ'^. 

A random vector field is i.i.d. iff these random variables are independent and iden- 
tically distributed. The components of a random vector field in a particular direction 
are a random field in the sense of Definition 12.1.21 ^ f^ict which we will use frequently. 
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We will say that a random vector field satisfies the Lyapunov condition if each of its 
components does in the sense of Definition 12. 1.31 

A random vector field fj is isotropically distributed iff for each x G Z"' and R G 
SO{N) the distribution of fjx is the same as the distribution of Rfjx- Among other 
things, this implies that the component e ■ fjx in an arbitrary direction will have an 
absolutely continuous distribution so long as fj ^ with probability one, and will have 
no isolated point masses (see the statement of Proposition 12.3.11 below) provided that 
ff ^ with nonzero probability. 

We then define systems by the quenched local Hamiltonians 

^hZ^ = ^ BriMX)) + Y.if' + Cx) • Br{K^), (2.1.22) 
X xer 

where each 

K/x ^ 

vector operator, that is a collection of N operators satisfying 

RK:c = Rx^KxRx: (2.1.23) 

and we also assume that the components of Kx are in Ax- Other local Hamiltonians, 
free energies, etc. are defined in the same terms. Then if ff is isotropically distributed 
and ^'o is invariant, we will say that the system described by H^'^ is isotropic. 

We will have need of a restriction on long range interactions similar to Assump- 
tion [2TLS] to extract additional results for these systems. The assumption (employed in 
the proof of Lemma l4.5.ip is as follows: 

Assumption 2.1.11. The sum 

^(diamA:)2|A:|||^'o(A:)|| (2.1.24) 

is finite. 

For pair interactions, this reduces to the statement 

5^ ||M/o({0,x})||x||^ <oo (2.1.25) 

x€Z^ 

found (in slightly different notation) in 
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2.2 Thermodynamic limit and notions of long range order 

The first requirement in talking rigorously about the thermodynamics of an infinite 
system is to prove the existence of the thermodynamic limit of some basic quantity. 
For lattice systems one conventionally uses either the pressure (as in [471 |46j ) or the 
free energy density (as in [S^iSGj) - they are related by P = — /?/, so for most purposes 
they are interchangeable. We will employ the free energy density, since it has the 
considerable advantage of having a well-defined behavior at /3 = oo (zero temperature) 
where in the absence of residual entropy it coincides with the ground state energy 
density. 

We define the free energy density for a finite system in the more or less obvious 
manner, as 

F/?r(C,w) 

ftBiC,uD:= , (^-^-^^ 

where |r| is the number of points in T. As the notation suggests, this depends on the 
choice of boundary conditions and of the disorder variables. In the thermodynamic 
limit, however, the dependence on boundary conditions disappears and the dependence 
on the disorder variables becomes trivial, as the following theorem will show. Essentially 
the same statement was first proven by Vuillermot in 1977 [62j : the version given here 
is more suited to the present work. 

Theorem 2.2.1 ([3^). Let Assumption \2. 1 ."U be satisfied. For any h, any i.i.d. random 
fields ri,v with finite variance, any (3 G [0,cxd], there is a set N G S^^"^ such that 
^ [{'Hi y) & M] = I so that the limit 

^(/3,/i) := hm /r\,B(C,a;) (2.2.2) 

exists for all {C,,iS) G N , /i G and all B, and is independent of Q, w, and B. 
Furthermore, 

lim " ^ , " = (2.2.3) 

for all (C,w) G M. 

This theorem was stated for classical systems, but the proof depends only on some 
properties of / - in particular Lemma l2.1.5l - which also hold for quantum systems. 
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The fact that the hmiting free energy is almost certainly independent of the random 
field provides the following conclusion: 

Corollary 2.2.2 (Brout's prescription [2l [62]). 

hm Avf^^^siCui) = Av hm /r\,B(C,a^). (2-2.4) 

In other words, one can take the average over the randomness before or after the 
thermodynamic limit without changing the free energy. 

Since is a limit of convex functions, the following useful fact (also noted in [3]) 
follows immediately from Theorem 12. 2. II 

Corollary 2.2.3. T{f3,h) is convex as a function of (3 and concave as a function of h. 

This allows us to prove some handy results which extend the relationship between 
the derivatives of the free energy to expectation values of certain observables from finite 
to infinite systems. To begin with, note that 

-Qh " M ^ ^'^'■^)t,b (C, ^)- (2.2.5) 

The convexity of J- does not imply that the above derivative always converges in the 
thermodynamic limit, but it does imply something almost as good: 



Corollary 2.2.4. 



L— >-oo Fr, 



dh- ' dh+ 



(2.2.6) 



where LIM denotes the set of accumulation points, and denote directional deriva- 
tives with respect to h. 

The above statement is about the average of (k) over the whole system, or in other 
words it is a statement about "long long range order". It is also possible to make a 
similar statement relating to "short long range order" : 

Theorem 2.2.5. 

dF dF 



LIM LIM i^^)kuBi(^^)^ 



dh- ' dh+ 



(2.2.7) 
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Proof. Let F^'^^ denote the free energy with the fixed field within A changed by 6, so 



that 



1 V- . 1 dF!i'^'% 



(2.2.. 



<5=0 



Now from Lemma 12.1.51 we see that 



' ^^VJi - FVIi) = /r^o - fLo + O f ^ ) (2.2.9) 



uniformly in M. Then for (C, € J\f, this implies that 

1 / ^h,5,rL TphATL 



^ rrM,i^ - ^rTB" j = h + 5)- ^(/3, /.) (2.2.10) 
and the conclusion follows by standard convexity arguments. □ □ 

Choosing a positive sequence (5i — )• such that F is differentiable at all h ± 5i, we 
have also 

dF 



lim lim lim — - V {k^)'^'^% {( , uf) = lim — 

i^-oo L^oo M-i-cxD 1 r ^ — ' L,^ z^oo Ofl 



h±s, dh± ' 



(2.2.11) 



which together with the individual ergodic theorem (applicable since the random fields 
are i.i.d) this implies 

Corollary 2.2.6. 

B F 

lim lim Av(«,.)^fj5 (C,^) = J—, (2.2.12) 

2.3 Statement of main results 

The first main result of the following chapters is: 

Proposition 2.3.1. In dimensions d < 2, any system of the type described in in 
Section \2.1{ with r] an i.i.d. Lyapunov random field and 7 i.i.d, has F differentiable in 
h for all h, provided any of the following hold: 

• The system satisfies the weak FKG property with respect to k, j3 < 00, and the 
distribution of r] is nontrivial 

• (3 < 00, and the distribution of 7]q has no isolated point masses (i.e. there are no 
real numbers x and 5 > Q such that F [|??o — x\ < S\ = ¥ [rj^ = x] > 
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• The distribution ofr]Q is absolutely continuous with respect to the Lebesgue measure 

We note that the phrasing of the result relates to the way we have arranged the 
Hamiltonians of the systems under consideration, so that the result has something to 
say only when the random field can be expressed as part of the source field for the 
order parameter, in other words when the randomness is direct in the sense used on 
p. [5] above. 

We also establish 

Proposition 2.3.2. In dimensions d < 4, any isotropic system of the type described 
in Section \2.1.1\ satisfying Assumption \2.1'AT\ has V^J-" continuous at 0, provided the 
distribution of fj is isotropic and one of the following holds: 

• |r/o| > with probability 1, or 

• (3 < oo, and the distribution of fj is not concentrated at a single point. 

We note that an apparently weaker condition on the distribution of fJ is adequate 
because the what will ultimately be important is the distribution of a particular compo- 
nent. With the assumption of an isotropic distribution for the vector, the components 
satisfy the stronger conditions used in Proposition 12.3. 1] or Theorem 13.3.21 as discussed 
above. 
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Chapter 3 
A nonlinear central limit theorem 

3.1 Background 

The "classical" central limit theorem, long one of the central elements of probability the- 
ory, states that a sum of independent random variables with finite variance converges, 
in distribution and on an appropriate rescaling, to a normally distributed random vari- 
able. There are a number of generalizations, the best-known due to Lindeberg |63j . 
which generalize this notion by replacing the i.i.d assumption with a weaker assump- 
tion, including the possibility that the distribution of the variables, as well as their 
cardinality, changes as the limit is taken. 

We will present a clarified version of a result due to Aizenman and Wehr [3j which 
builds on results of that kind to replace the customary sum with a member of a much 
larger class of functions, which however have certain properties (a partial symmetry with 
respect to permutations of arguments, and a fairly strong continuity) in common with 
it. This exposition also incorporates a necessary correction pointed out by Bovier [56]. 

We should note that the statement that a certain sequence of random variables, 
described as a related collection of Lipschitz continuous functions of a family of in- 
dependent random variables, converges in distribution to a normal random variable, 
is closely related to the concentration of measure phenomenon [64] . Among its many 
other facets, this involves upper bounds on the probability with which certain classes 
of random variables described as functions of a family of N independent random vari- 
ables. A central limit theorem involves an estimate of a similar form. The result we 
will discuss is stronger than a concentration estimate in that it provides a lower bound 
as well as an upper bound (the latter would not be useful for our main result); however 
it is purely asymptotic, whereas concentration of measure techniques provide speed of 
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convergence information as well. The assumptions are in some ways stronger and in 
some ways weaker than those involved in concentration of measure: 

1. The result below uses a form of Lipschitz continuity based on the £i norm, which 
is stronger than the £2 notion used in concentration of measure (see below) but 
more suited to functions of infinitely many variables. 

2. The result below assumes translation covariance (a weak form of exchangability) , 
but no assumption is made on its level sets. Gaussian random variables do not 
play a distinguished role. 

3.2 Definitions 

As well as the notions related to random fields introduced in the previous chapter, we 
will make use of the following: 

Definition 3.2.1. A function r : Z'^ x £^ — >■ R (equivalently, a collection of functions 
of fields indexed by elements ofU^) is translation covariant ifTx{rf) = Tx-y{Tyr)) for all 
x,y,ri, where Ty denotes translation. 

We will have occasion to frequently use the norm on S, 

llClli := E IC.I; (3.2.1) 

in particular this defines a Lipschitz seminorm on functions / : f — > M by 

Ill/Ill := sup (3.2.2) 
CC'eS IK - 4 111 

0<||C-C'lli<oo 

It is worth spending a moment on the comparison of this norm with the similar quan- 
tity based on the i'^ norm which appears frequently in the concentration of measure 
literature. The £^ norm in this context is defined by 

IICIl2:=fEC') , (3.2.3) 

and the related Lipschitz seminorm by 

111/1112:= sup (3.2.4) 
CC'G^ lis - S II2 

0<|IC-C'll2<OO 
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In concentration of measure one is concerned with functions of N variables, for which 
the supremum in the above expression is attained with and C' differing only in the 
corresponding N elements, whence 

||C-Clli<^||C-C'||2 (3.2.5) 
by Young's inequality. Using this to compare Equations (|3.2.2p and (j3.2.4p . we see that 

|||/7v||l2<^lll/7v|||, (3.2.6) 
so that if we have |||/Ar||| < 1 (as below), this implies |||/Ar|||2 < VN- 

3.3 The proposition 

Proposition 3.3.1. Let rj be an independent, Lyapunov random field, and let Gl '■ 
f — )■ M 6e a family of functions indexed by L £ N, each with the following properties: 

1. Gl depends only on the values of the field for sites inTi 

^. IIIClIII < 1 

3. AvGiir]) = 

4. Av [GLivMriA = Ca] = Gl' o whenever T^x^u = A C 
Then 

GUr])/L^/^ ^N{0,b^) (3.3.1) 
in distribution as L ^ 00, for some b satisfying 

AvG? < 6^ < 2(Av|7?o|)^ (3.3.2) 

In order to use this result, we will need to establish some control over the conditions 
under which Av > 0. To do this, we employ the following theorem (proven in 
Appendix III of [3]): 

Theorem 3.3.2. Let v be a Borel probability measure on M, and 

{geC\W)\\\\g\\\<l,\\\g'\\\<fi}, /3 < 00 
■= \ , [6.6.6) 

{geG{R)\\\\g\\\<l}, /3 = 00 
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and also 




n 1/2 . 

g{xfu{dx) g e Vi,;3, / g'{x)iy{dx) = M 
g(^Vi,p,g'{-)>0,J g'{x)v{dx) = 



.1/2 



g{x) u{dx) 



M 



(3.3.4) 
(3.3.5) 



1. e,(0,/3) =7,(0,/3) = 

2. 9u{M,0) is nonzero for all M > iff i' is absolutely continuous with respect to 
the Lebesgue measure 

3. For finite {3, 9,^{M,f3) is nonzero for all M > {) iff u has no isolated point masses 

4- For finite j3, 'j^(M,P) is nonzero for all M > iff is not concentrated at a 
single point 



To employ this, we note that Gi £ Vi^b for 

dGi 



B :-- 



dCo 



(3.3.6) 



if the derivative on the right hand side exists everywhere, and B = oo otherwise. Then 
when Gi has a distributional derivative G'^ with Av G'l = M for some M > C0 we have 
AvG? > 0l{M,B) (for monotone Gi, AvGf > j'^{M,B)). 

We will not provide a proof of Theorem I3.3.2t but since it is the source of a per- 
plexing limitation in our result (as in the classical case) some commentary seems to be 
warranted, and it is possible to provide some insight into the situation and its prospects. 
This will be done in Section 13.51 below. 

Before going on to the proof, we should clarify the relationship to the formulation 
of the corresponding result. Proposition 6.1 of [3J. Much of the difference is due to 



^Note that G'l is a function of only one variable. The existence of a distributional derivative, i.e. a 
Lebesgue-integrable function satisfying Gi{x)dx — Gi{y) — Gi{x) is guaranteed by item[2] which also 
implies that |[G'i||oo < 1 and therefore also that AvG'i exists. It may be, however, that this derivative 
is not unique, and when the distribution of 770 is not absolutely continuous with respect to the Lebesgue 
measure and B = cxd it is possible that this could allow more than one valid choice of M, although this 
is immaterial for the application we have in mind. 
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the fact that I have separated the main result from the positivity criteria embodied 
in Theorem 13.3.21 but there is a remaining difference in language in which the results 
are framed, as the following lemma, which also brings the abstract objects of Propo- 
sition [3]3TT] into a form more closely related to their use in a thermodynamic context, 
should clarify: 

Lemma 3.3.3. Let ij be an independent Lyapunov random field, and let Gl :<?—)■ M 6e 
a family of functions indexed fey L € N, and : — s- M a translation covariant family 
of functions satisfying 



2. AvTx = M for some M > 



0, x^Fl 



(3.3.7) 



3. \TxiC)\ < 1 and 
I AvGL(r/) =0 



< B' for all X, Q 



Then r] and Gl satisfy the hypotheses of Proposition \3.3.I\ (with the same M). 
Furthermore: 

1. If Tx is nonnegative, then Gi is nondecreasing. 

2. Gl G Vi,B' 

3. Gl has a distributional derivative G'l with Av G'l = M 

Proof. The first three conditions in Proposition l3.3.l] are trivially satisfied since 1 1 IG^^I 1 1 < 
sup^ kxlOl ^ and G'l = Av [To|ryo]. The enumerated properties of Gi are equally trivial. 

For the last point, we note that Av [GL|?yA] and Gl' o have derivatives given by 
identical expressions in terms of r^,., and so can only differ by a constant; but both have 
zero mean, and so that constant must be zero. □ 



33 



3.4 Proof of Proposition 13.3.11 

Let Gl and r] satisfy the conditions of Proposition 13.3.11 Order the elements of JJ^ 
lexicographically, and let F{L, k) be the set consisting of the first k elements of F^, (of 
course < A; < L'^). Then we define 

Y^k ■= Av [GL\r]F{L,k) = CFiL,k)] - Av [GL\vF{L,k-l) = CF(L,k-l)] (3-4.1) 

so that 

Gl{C) = Y,YlMC) (3.4.2) 

k=l 

for all C G £■ The definition l3.4.1l of Y makes it a martingale array and we will ultimately 
obtain a proof by showing that it satisfies the conditions of an existing central limit 
theorem for such objects j65J . As is the case with other central limit theorems for 
non-i.i.d. arrays, the conditions of this theorem are basically the existence of a limit of 
the average variance (Lemma 13.4.41 below) and the vanishing of fluctuations on a larger 
scale (Lemma 13.4.51 below) . 

The following result tells us that fluctuations in Yi^k{r]) (a random variable) are 
basically no worse than those of the field at a single site. From here on, we will let 
Xk denote the fcth element in T^,, when the value of L is clear from the context, and 
Cfc = Cxk etc. 

Lemma 3.4.1. For all C, & £, 

|^L,fc(C)l < (|a|+Av|r?o|). (3.4.3) 

Proof. We can write 

^L,fc(C) = Av (^GL(r/pc,Ci,...,fc,%+i,...,Ld) - G'L(7?rc,Ci,...,fc,%+i,...,2.d)) ; (3.4.4) 

then the assumption that |||Gi^||| < 1 means that the quantity being averaged above 
has absolute value no more than | — "Hk]- Thus 

I^L,fc(C)l < Av ICfc - ilk\ < (ICfcl + Av |r?o|). (3.4.5) 

□ 
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This is a uniform bound in absolute value by a square-integrable function (since the 
Lyapunov condition implies in particular that r]x has finite variance), and will allow 
us to apply a number of general convergence theorems. In particular it makes any 
collection of the functions Yj^^k uniformly integrable, and we will take advantage of 
this to show that the asymptotics of Y are described by a translation covariant (from 
another perspective, stationary or exchangeable) object W. 

For X G Tl, let Yl,^ denote Y^^k for k such that = x, and let = J^l,L''- Then 
it is evident from the consistency condition on G^. in Proposition 13.3.11 that 



Yl^x = Av [Yl'^xIvTl = Ctl] > 



(3.4.6) 



which is to say that for fixed x, the sequence Y^^xiil) forms a martingale with respect 
to J-, and applying the uniformly integrable martingale convergence theorem we have 

Corollary 3.4.2. For each x G Z'^, the C} limit Wx = liraL^ca Y^^x exists with 

YL,x = AY[Wx\vr^=Crr.] (3-4.7) 



whenever B x. 



Lemma 3.4.3. Wx form a translation- covariant family. 



Proof. Recalling the definition of translation covariance (j3.2.ip , we examine 



Wx{TyC) =C^-\imYL,x{TyO 

L-^oo 



: £^-lim ( Av 



G 



L+\\y\\ 



•T.„ 



'nT-yF{L,k) — CT^yF{L,k) 



(3.4.8) 



- Av I^G^^+llyll^ O Ty VT-yF{L,k-l) " F(L, fc - 1 ) 

where we have obtained the right-hand side by writing out the definition of Y^^k and 
some elementary properties of the conditional expectation. We then rewrite the right 
hand side again using the consistency assumption on Gl , Assumption 2] of Proposi- 
tion [3i3?Tl to change coordinates, and obtain 

Wx{TyC) =£/-limAv Y, 



L+\\y\\oo'^-y 



VT-yFL = Cr-yTL 



CUimAv [Wx-y\vT.yr, = Cr.^rJ = Wx-y{0 



(3.4.9) 



□ 
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We are now ready to prove 



Lemma 3.4.4. Let 6^ = AvW^; then 



1 ^ 

^ Av [Yl,,\r]p^L,k~i) = CF(L,k~i)] - 



k=l 



(3.4.10) 



(3.4.11) 



in measure (and therefore also in distribution) as L — t- cx). 

Proof. We will do this essentially by showing that Y^^x can be replaced by Wx, apart 
from a boundary term which vanishes in the limit L — t- oo. We can of course write the 
summand above as 

Av [Ylk\r]F{L,k~i) = CF{L,k-i)] = Av [W^^\r]<xk = Ux^] 

+ Av [W^^ \r]F{L,k-l) = CF(L,k~l)] 

-Av [Wl\v<x, = C<x,] 

+ Av [Ylf, - Wx^\r]F(L,k-l) = CF{L,k-l)] , 
where by ?/<Xfc and similar expressions we mean rj^ for x < in the lexicographic order; 
we then deal with the different terms separately. Letting /(r?) = Av [VFq |7/<o = C<o] i we 
use translation covariance and the fact that the conditional expectation is a projection 
in to obtain 

\\^y[Wl\vFiL,k-i) = CF(L,k-i)] -Av[Wl\v<x,]\\2 ^3^^2) 
<||/-Av [/|r?r, = CrJ|l2=:«i(^) 
where R is the largest integer for which Txf.Tji C T^. Employing Holder's inequality 
followed by a similar step, we have 

II Av [Ylf, - W^Jr]p(^L,k-i) = CF(L,fc-i)] 111 

< \\YL,k + WxJ^ II Av [YL,k - Wx,\7ip^L,k-i) = CF(L,fe-i)] II2 (3-4-13) 

< 2 \\Wx, II2 11^0 - Av [Wo\r]rn = CrJIls =: a2{R). 

The i2^-norm expressions used to define oi and 02 must vanish as i? — )• 00 and depend 
on L only through R since / and Wq are square-integrable. Since the relevant terms 
in Equation (13.4. lOp are an average over k in which the proportion of the terms with 
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arbitrary large R increases without bound as L increases, these terms go to zero in 
measure. The proof will be complete if we can show that 



1 ^ 

^ Av [W,, I r,<,, = C<. J ^ Av 



(3.4.14) 



k=l 



in measure; which, given the translation covariance of Wx and the fact that r] is i.i.d., 
follows immediately from the £^-ergodic theorem. □ 



To obtain inequality ()3.3.2p . we note that Lemma 13.4.11 implies a similar bound on 
\Wq\, and therefore that AvH^q < 2 Av |77oP, and that 



Av > Av ( Av [Wo\r]o = Co 



(3.4.15) 



By dominated convergence of conditional expectations (applicable by Lemma [3.4.ip and 
the definition of Y in Equation ()3.4.ip , 



Av [Wo\r]o = Co] = lim Av [Yl,o\vo = Co] = Gi, 



(3.4.16) 



and so 

AvG? < AvWo^ < 2(Av|7?o|)^. (3.4.17) 

All that remains is to show that we have a sufficiently strong control on the large 
fluctuations of YL^kiv)- 



Lemma 3.4.5. For any a > 0, 



VF{L,k-l) — CF(L,fc-l) 



(3.4.18) 



fc=l 



in probability as L oo. 



The proof is due to Bovier [56], and provides a correction of a mistake in [3]. 

Proof. Note that the average of the left hand side above is 
1 



Av 



1 

1? 



k=l 

^Av fel > 



\YL,k\ > aV^I"^ VF{L,k-l) = Cf(L,A:-1) 



1 



k=l 



(3.4.19) 



r2q 

'L,k 



1/g 



P 



\YL,k\ > aL'^/^ 



i/p 



k=l 
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for any \/p + l/q = 1 by Holder's inequality. Chebyshev's inequality (noting KvY^^k 
0) gives 



P 



< 



L,k 



(3.4.20) 



and Lemma 13.4.1 1 together with the fact that rj^ has finite variance gives a uniform (in 
k) upper bound on AvY^^, so the right hand side goes to zero; we could conclude 
that the right hand side of (13.4. 19p is zero if Av Y^'^ is finite for some finite p, i.e. for 
q > 1. Lemma 13.4.11 implies that AvY^^ is finite if rjx^ has a finite 2q moment, and 
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the Lyapunov condition is precisely the fact that this is true for some (7 > Icl We then 
have 



lim Av 



1 

I? 



k=l 



\Yrk\ > aL'^'^ 



'>]F{L,k-l) — CF{L,k-l) 



0, (3.4.21) 



which, since the quantities inside the average are a uniformly integrable family of non- 
negative functions, can only be true if those functions converge in probability to 0. □ 



We can now apply Theorem 3.2 of [65] and conclude that 



GL(r?)/L'^/2^iV(0,62), 



and the proof of Proposition 13.3.11 is complete. 



3.5 Bounds on 6^ 



(3.4.22) 



Theorem l3.3.2l is the source of a perplexing limitation remaining in our result, the exclu- 
sion of distributions with isolated point masses (in particular of discrete distributions) 
from most of the result. Since the proof of this theorem is rather opaque, it seems 
worthwhile to give a heuristic discussion which may make the result appear less arbi- 
trary, and clarify some of the issues involved in attempting to obtain a more powerful 
result. 

Consider the simplest nontrivial discrete measure: let v = + Given M and 

/?, is there a function g satisfying the apparently relevant properties of Gi (111(7111 < 1, 



■^This is the only place where our results require the full Lyapunov condition, and not merely existence 
of 2 moments. 



38 



fg'di^ = M, \\\g'\\\ < /3) with f g^di^ = 0? 
whenever M < mm(/3/3, 1), the function 



Quite often the answer is yes. For example 



9i{x) = < 




1) 



x), -l<x<l 



X > I 



(3.5.1) 



M{x + 1) 



X < -1 



clearly fits these requirements. This is of course not a monotone function, and indeed 
it is more or less obvious that (as Theorem 13.3.21 states) no monotone, continuously 
differentiable function will do: if ^'(±1) > 0, then for such a function g{^l) ^ 0. 

This would appear to leave little room for improvement in Proposition 13.3. H but 
this is not quite the case. Proposition 13.3.1] stipulates that Gl{C,) should be monotone 
in C,o for all L and but all that is needed is that Gi be monotone, which should 
be a weaker requirement. It is, however, not immediately clear that it follows from 
something as well-known as the FKG inequalities which imply the special case. 
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Chapter 4 
Free energy fluctuations 



4.1 Definition of Gl 

In order to prove Proposition 12.3.1] we need to construct a sequence of functions Gl 
which represent the effect of the random field r/ on the free energy difference between 
states with the largest and smallest permissible values of the order parameter. Since 
our most robust way of accessing these states is by taking a limit in the uniform field h 
which couples to and since we are concerned with the thermodynamic limit, it should 
be plausible that one candidate is described by the formal expression 

1 



lim lim — Av 

5-J.0+ Af->oo 2 



(4.1.1) 

where ri is the function which sets the field to inside Tl; Fm is the free energy 
of the system on the finite domain with periodic boundary conditions (defined in 
Equation (|2.1.5p ). and for brevity we have omitted the argument v, or in other words 
we let F(C) = Av [F{ri,v)\r] = and rjr^ is the collection r]x where T^Aq C Tl and 
likewise for similar expressions. 

However it is hardly clear that the expression (14.1. ID is well-defined. We will show 
that a quite similar quantity is, but first let us turn to a few observations which should 
help motivate this choice. 

For convenience, we let 



GiMO ■■= \ Av [F^t\ri) - Fjl+'iTLiv)) - Fj^-'iv) + Fj^-'irUv)) 



(4.1.2) 
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for M > L. Then 



Av 
Av 



h-S 
M 



drjx drj. 



(4.1.3) 



This means that 



r"5 



< 1 uniformly in all parameters, which will carry over in 



the limit M — t- oo to assumption [2] in Proposition 13.3.11 It also means that 



Av 



which in light of Corollary 12.2.41 should mean 



i Av ((/€x.)m iv) - {kx)m ^ iv)) 



(4.1.4) 



(4.1.5) 



2 \dh+ dh- 

allowing the desired control on h. The use of a conditional expectation in the defini- 
tion should take care of assumptions [U and [U and we can arrange for the remaining 
assumption (mean zero) by simply subtracting the mean value. 
Let us return to this more carefully: 

Proposition 4.1.1. Let at least one of the following hold: 

1. p <oo 

2. The distribution of rjx is absolutely continuous with respect to the Lebesgue mea- 
sure 

Then there is a decreasing sequence (5j — >• and an increasing sequence of integers 
Mj — >• oo such that 



Gl{C) := lim lim (g'Im -Av 



L,h'U 



(4.1.6) 



exists for all Q ^ £ and all L ^'H. Furthermore, the family Gl satisfies the hypotheses 
of Proposition \3. 3.1\ and Gi has a distributional derivative G'l satisfying 

I f dF dF 



AvG'i 



(4.1.7) 



2 \dh+ dh- 

I wish to point out that the proof will not assume that r] and v are mutually 
independent, but only that the different t]^ remain independent when conditioned on v 
- this will be important for systems with continuous symmetries, where t]^ and will 
represent different components of a random vector, and may therefore be correlated. 
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4.2 Proof of proposition l4.l!T] - finite temperature 

In this situation we will proceed by constructing functions r which satisfy the conditions 
of Theorem I3.3.3[ 

The following convergence argument will be used frequently in what follows. 

Lemma 4.2.1. Let fij : M.'^ — )• M 6e a family of functions labeled by i,j G N, each 
satisfying < 1 and /ij(0) = 0. Then there are subsequences ik,ji such that 

/(z)= hm lim (z) (4.2.1) 

exists for all z G M^. Furthermore the convergence is uniform on any compact H C M.'^ . 

Proof. Note that the condition < 1 implies uniform equicontinuity. On the 

compact domain H„ := [—n,n]^ we have the uniform bound \fxy\ < Nn, and so by 
the Arzela-Ascoli theorem, any infinite collection of these functions has a subsequence 
which converges uniformly on 

We then apply the diagonal subsequence trick as follows: there is a sequence jf so 
that fi ji converges uniformly on Hi, which has a subsequence jf so that /2 j2 and f^ji 
converge uniformly on H2 and so on. Then the diagonal subsequence ji = j\ has the 
property that for any k, n converges uniformly on i^ri, with limits : M — y M 
with the same properties we have used above. By the same argument, we can now 
choose a sequence ij, so that /ji converges uniformly on Hi, a subsequence i| so that 
/j2 converges on E2, etc. Then with i^ := i^, ji := j', we have the desired result. □ 

The same argument also gives 

Lemma 4.2.2. Let fij : — t- M 6e a family of functions labeled by i,j G N, each 
satisfying < 1 and \fij{z)\ < c < 00 for all z G R^. Then there are subsequences 

ik,ji such that 

/(z)= hm lim/,,,,(z) (4.2.2) 

fe— >-oo 1^00 

exists for all z G . Furthermore the convergence is uniform on any compact E C . 
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Now let 



niAO ■■= = I (Av iC,v) - Av iC,v)) ; 



(4.2.3) 



note that the quantity defined does not depend on L, and that C is fixed in the averages 
which are taken over the other fields v. To simplify the similar expressions appearing 
below we will write (•)^_^ {() := Av (•)^^ iC^ll)- Also, evidently 



?Xo(C)l < ll'^.ll = 1 



(4.2.4) 



and 









2 



{l^xl^y) ]^ (C) {KxKy)j^,^ (C) 



{^x)if (c) (0 + {^x)i7' (0 («.);;7' (o 



<2/3; 



then 



'Ai.A/.^lO := Av 
obeys the same bounds, 









drjx 



Av [OliAv)\vr, 



(4.2.5) 



(4.2.6) 



\^ImAO\ < 1 

< 2/3 



5Cy 



(4.2.7) 
(4.2.8) 



For f3 < oo, this means that for each L, cj)^ g is a uniformly equicontinuous family 
of functions of the L'^ variables ■ We can apply Lemma 14.2.21 to find a decreas- 
ing sequence 5j — )• and an increasing sequence Mj — )• oo (by applying the diagonal 
subsequence trick, we can choose them to be independent of x and L) so that 



rdC) ■■= lim lim </>i,M„5.(C) 



exists, and by uniformity of convergence 



GUC) := lim hm {G%,, {v)-AyGl,, 



also exists with 



OGl 
dCx 



(4.2.9) 



(4.2.10) 



(4.2.11) 



for all L, C, X & Tf^. 
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Note that Equation (j4.2.6p implies that, for any K < L < M , 



Av [ct>l 



L,M,5\'nrt 



(4.2.12) 



and by the conditional form of the dominated convergence theorem this implies that 



= Av [iI^Hvk], 



(4.2.13) 



which makes ^ martingale; < 1 makes it a uniformly integrable one, and 

applying the relevant martingale convergence theorem we see that 



exists as an limit, with 



:= hm Vl 



Av[r,,|7?rJ = V'!-^^^ 



(4.2.14) 



dCx 



(4.2.15) 



< 2/3. 



Following through the various limits, we see that |t^(C)| ^ 1 and 

We have defined in terms of periodic boundary conditions, so Equation (j4.2.12p 
also implies 



(4.2.16) 



which is translation covariance. Following through the limits used to define -0 and r 
(thanks to the fact that the sequences involved are independent of L and x), we see 
that this implies 

T,+y{TyC) = r,(C). (4.2.17) 

Finally, 



AvTa; = Avtq = Av0? = lim lim Av0? »r, = lim lim Av 



,0 



(4.2.18) 



and applying Theoroms 12.2.41 and 12.2.61 we have 



1 



Avr^. = lim lim ^ Av ( ^ (^) " ('^^)m ^ 



2 \dh+ dh- 



(4.2.19) 
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4.3 Proof of proposition l4.l!T] for absolutely continuous distributions 
of r] 

First of all, note that it is obvious from Equation 14.1.21 that G'^^(O) = 0; Inequal- 
ity ()4.2.4p implies < 1, so we can apply Lemma [4.2.1l to obtain (j4.2.10p . with 
1 1 Gill < 1. AvGl = is obvious. We can apply the diagonal subsequence trick to 
obtain sequences independent of L, which implies that the consistency condition H] of 
Proposition 13.3.11 is satisfied. 

The hard part is to show AvG'i = M. Without uniform equicontinuity of the 
derivatives, we have no reason to expect that an object like (j) of the previous section 
will converge uniformly, and without that we have no reason to expect that a pointwise 
limit will still be a derivative. However the following theorem allows us to find a 
particular kind of weak limit which will do the trick: 

Theorem 4.3.1. Let g be a measurable function and a sequence of measurable 
functions such that \\gn\\ < 1; ll^ll ^ 1; CLnd 

r-b fb 



lim / gn{x)dx = / g{x)dx (4.3.1) 

J a J a 



for all a,b € M. Then for any signed measure v with finite total variation which is 
absolutely continuous with respect to the Lebesgue measure, 

lim / gndu = / gdi^. (4.3.2) 

n-^ooj J 

Proof. Without loss of generality we can consider only positive measures (which are 
all that is necessary for the present work anyway), thanks to Hahn's decomposition 
theorem |6B]. As a preliminary, we see that for any Borel set A contained in a bounded 
interval I 



n— >oo 



lim / gn{x)dx = I g{x)dx (4.3.3) 
" 'a J a 



since for any e there is a set which is a finite union of intervals which approximates 
A in the sense that \{E^/\A) < e. Then also 

gndX- [ gndX <e (4.3.4) 

J A 
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uniformly in n, so we can take the limit e — )• and exchange the order of the limits to 
obtain Equation (j4.3.3p . 

There is a nondecreasing sequence am of simple functions so that am ^ pointwise, 
and each am has bounded support. By Equation ()4.3.3p . 

lim / gn{x)amdx = / g{x)amdx (4.3.5) 

n^ooj J 

for all m. In fact since am are a nondecreasing sequence, we can apply the Beppo Levi 
theorem [551 to obtain 

ml^oo J f« ^ / ^^^^"^ ^ / ^" (4.3.6) 



and thus 
Furthermore, 

dX - I gndu 



lim / gnCTmdX = / g-ndu. (4.3.7) 

m^ooj J 









~ dX, 






[dX 





— - CTm 1 dA 



(4.3.8) 



and since this last bound is independent of n, the convergence in Equation (|4.3.7p is 
uniform in n. Taking the limit n — ?• oo and exchanging the limits on the left hand side 
gives Equation ()4.3.2p . □ 

The set of absolutely continuous finite signed measures is isomorphic to £^(R), the 
predual of C°° (M) , so the substance of Equation (|4.3.2p is also expressed by saying that 

— ^ 5 in the weak-* topology of £°°(M). This is a convenient way of phrasing the 
following 

Corollary 4.3.2. Let fn be a sequence of functions M — )• M such that fn^f pointwise, 
o-'i^d |||/n||| ^ 1- Then their distributional derivatives converge to the distributional 
derivative of f (f^ — t- f) in the weak-* topology o/£°°(M). 



comparable statement appears in the proof of Rademacher's theorem in thus the proof here 
is more logicaUy circuitous than necessary, but we hope it is the most intelhgible way to convey things 
to our readers. 
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Proof. Thanks to Rademacher's theorem [67], Lipschitz continuity guarantees that the 
distributional derivatives, i.e. functions satisfying 

rb 



f'^dx = fn{b) - fn{a) (4.3.9) 

J a 

for any a, 6 E M, exist with ||/^|| = |||/n||| < 1, and the convergence fn^f then imphes 

fb rb 

hm / f'r,dx= / f'dx. (4.3.10) 

J a J a 

This allows us to apply Theorem 14.3. II and the result follows immediately. □ 



Applying Corollary 14.3.21 twice to G"J'^j.(r/) — AvG^'^j, and Gi, we obtain 



AvG'i = lim lim Av {G\' 



(rj) 



(4.3.11) 



and by Equation (j4.2.3p and Theorems 12.2.41 and 12.2.61 we obtain 
and the proof of Proposition 14.1.11 is complete. 



(4.3.12) 



4.4 Proof of Proposition 12.3.11 

We now turn to the boundary estimate ()2.0.ip . Let A^ be the smallest subset of 'E'^ so 
that TxAq C Al (i.e. Kx S •Aal) for all x G Tl, and for Tm ^ A^ let 



^m|l(C) - Av logTr exp [-PH'/J - 



(4.4.1) 



where the subscript refers to free boundary conditions, and the subscript * refers to 
periodic boundary conditions on the edge of rj\/ and free boundary conditions on the 
edge of A^; this lets us write 



h,C,,uj_ Trh-Xtiil 



Th,C,,ui 



f\rL*' 



(4.4.2) 



whence, by Lemma 12.1.51 



< Av 



V, 



(4.4.3) 
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then 



■h~S/ 



?h-S/ 



< 



+ 4Av 



V, 



(4.4.4) 



Since H^^^ and H^'^^^^^ act on disjoint subsets of the lattice, they commute, and 



^M|L(C)=<,0(C)+i^AAL*(O. 



(4.4.5) 



where 



1 



Flf\L{0 := - Av-ilogTVexp ("KIaa,* 



(4.4.6) 

When we use this to expand the right hand side of (|4.4.4p . the Tm \^l terms cancel: 



F\ 



mli^) - Flt\li^rUri)) - F^f^lirj) + F-^(r^(r/)) 
and we can apply Lemma 12.1.51 again to bound this quantity, obtaining 



?h+5 



h-5/ 



(4.4.7) 



Futiri) - Ft^firUv)) " i^^^n '(^) + ^mT'(^l(^)) < 25|Al| = 0(61% (4.4 



M\L 



M\L^ 



where the last term is the effect of the constant field inside A^. 
Plugging this back into Inequality 14.4.41 gives 



Ft^Hv) - F^+'{rUv)) - F'-'iv) + i^i^r' (^L (^) ) 



which gives 

|G'i,M(C)l <2Av y; 

and since the 5 term is uniform in M, 

|Gi(C)| <2Av 



< 4Av 



v: 



+ 0{5L'^), (4.4.9) 



(4.4.10) 



V. 



(4.4.11) 



which we have assumed (Assumption I2.1TU|) to be 0{L'^ ^) + 0{L'^^'^). 

Now we need to show that this is in contradiction with Proposition 13.3.11 unless 
6 = 0. To demonstrate this absolutely clearly, we will convert these to statements 
about the moment generating functions of G^, Ave**^^. The distribution of a random 
variable is uniquely characterized by its moment generating function provided this is 
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finite on a sufficient region [66j , which in this context is guaranteed by ()4.4.1ip ; and then 
convergence in distribution is equivalent to pointwise convergence of moment generating 
functions. Then the conclusion of Proposition 13.3.1] can be restated 



lim Avexp ( tGi 

L—>oo 



exp{t^by2). 



(4.4.12) 



At the same time, if |Gl| < AL'^/'^ then 



(4.4.13) 



for all positive t; clearly ii b ^ 0, this will be incompatible with (j4.4.12p for sufficiently 
large t. Finally we note that Theorem 13.3.21 states that 6 = implies M = under any 
of the cases listed in Proposition 12.3. H and the proof is complete. 



4.5 Systems with continuous symmetry: Proof of Proposition 12.3.21 



As noted above, the main requirement of the proof of Proposition 12.3.2] is based on the 
improved bound 

\Gl{C)\<KL''-^ (4.5.1) 



which should hold at /i = We first note that Proposition 14.1.11 holds for the vector 
case with the following definitions, corresponding to Equations ()4.1.2p and (j4.1.6p : 



1 



GTMO = 2 FiUff) - Fll{rr.{fi)) - F^/^(r/) + F^/^(rz.(r/)) 



e-f]L = C,L 



G1(C)= lim lim (GiVrC)-Av 



G 



(e-ff) 



(4.5.2) 
(4.5.3) 



J— >oo J— >oo 

where e, an arbitrary unit vector, defines the component of the order parameter being 
examined. 

We can obtain the desired bound by focusing on 



5iV(C):=Av Fil{fi)-F^rifJ) 



since 



Gi(C) = -lim lim <M(C)-Av5tM(e-r?) , 



(4.5.4) 



(4.5.5) 



it is easy to turn uniform bounds on |5ffjvf(C)l into similar bounds on \Gl 
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Lemma 4.5.1. With g defined above, Assumption \2. 1 .U] imph 



les 



5iV(C)l=0(L^-^) 



d-2\ 



(4.5.6) 



Proof. Let p be the generator (in so{N)) of a rotation in a plane containing e, and 
for each x S Z'^ let px be the generator of the corresponding rotation in the single-site 
algebra We introduce the slowly varying angles 



0, 



L^ip^TT, < diix) < L 



(4.5.7) 



||x||i(ij;,(x) > L 

where diix) is the distance from x ioTi in the largest-component metric, i.e. 



diix) := min \\x - y\\^ . 
We also introduce the associated rotations on fields and on A defined by 



(4.5.8) 



/ X 

R9= (S) e^'"'- 

xeZ'i 



(4.5.9) 
(4.5.10) 
(4.5.11) 



Rg is unitary, and so we can rewrite the free energy Fj^j '^{fj) appearing in (j4.5.4p as 



Fm'{0 = --^logTrexp {-(iR^^ H-''^^^^Re 
we wish to use this to obtain something of the form 



Av 



Av 



logTYexp(-/3[/7j'^'^'^ + A/7e] 



e-f]L = CL 



which by Lemma 12.1.51 implies 



(4.5.12) 



(4.5.13) 



(4.5.14) 



^Unlike in 3], we will use the "mathematician's" convention that rotations are given by e*'', so p is 
an antisymmetric matrix and px is an antihermitian operator. 
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however this will not quite be sufficient, since we are not able to establish suitable control 
over ||Affe||. Instead, we will split F^^^^C,) in half and rewrite each part separately by 
applying an opposite rotation, to obtain 



Av 



e-ffL = CL 



1 

'2^ 



Av 



log Tr exp ( -I3[H^^'^'- + AHg 



+ log Tr exp ( + AH_ 



e-f]L = CL ■ 
(4.5.15) 

Combining the Cauchy-Schwarz inequality, the Golden-Thompson inequality, and 
Lemma 12.1.51 we quickly derive the general inequality 

Tre^_ \ 

(4.5.16) 



log Tr - log Tr e^^"^ - log Tr e*^/^ = log 



Tre^/2Tre^/2 



A 



B + C 



for arbitrary Hermitian matrices A, B, C. Applying this to Inequality (j4.5.15p gives 



CAf(C)< 2ll^^^ + ^^-^ll• 



(4.5.17) 



Now recah Equation ()2.1.22p : 



R.lH^'^^'^Re = Rg^ ( J]Pr(^o(^)) + J^^Se + Cx) • Pr{d^)+] Re (4.5.18) 
Since k are vector operators (recall Equation (|2.1.23p ). 

Cx ■ (Re^K..Re) = Cx ■ R9{^)x = [Re\C)x] ■ ^x\ (4.5.19) 

now inside there is no rotation, and outside we are performing an average with 
respect to an isotropic distribution, so this term makes no contribution to A//. 
As for the fixed field terms, we have 



e • (Rq ^KxRe) = {Rx ^e) • Kx- 



(4.5.20) 



The choices of p and 9 were intended precisely to make Rxe = —e for diix) > L; and 



for the remaining (3L)'^ sites we have 



e- (Rf) ^KxRe ] + e- Kx 



< 2, so these terms 



make a contribution to AH which is uniformly bounded in norm by 2(3L)'^5. 
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We are left with the terms arising from the transformation of the nonrandom inter- 
action. For any X and any (arbitrarily chosen) x G X n F, 



pzJdz~ex)pz 



(4.5.21) 



e-iey-e-)py Pr(^o(x)) (g) e^^^-'^^'^P^ , 
yyexnr J \zexnr J 

(using the rotation invariance of ^'o). Expanding the exponentials, we obtain 
R.ePT{^o{X))Re = Pr{^o{X)) + ^ (0:. - Oy) (pyPri^oiX)) - PT{^o{X))py) 

+ o((5l5=^ll*„(x,i| 



(4.5.22) 

where the estimate of the higher order terms uses 

„ , ttIIx — -hIIoo vrdiamX 
\0x - 0y\ < " ^ " < ^ (4.5.23) 

and the observation that the nth order term in the expansion is potentially a sum of 
terms, as well as ||Pr(^o(^))ll ^ ll^o(-^)ll- The first order terms are odd in 9, and 
will cancel in AHg + AH-q, with the leading term being second order. What appears 
there is 

Y: {R-oPrmx))Re - Pr(^o(X))) = O f ^ M ^'^''"f^^''"''' "''"^^^"l 
xnT^(!> \ xbo ' ' / 

= 0(L'^-2), 
(4.5.24) 

where the last equality invokes Assumption 12.1 .TTl 
Then we indeed have Equation (j4.5.15p . with 

\\AHe + AH_e\\ = 0{L'^-^) + 0{6L'^)). (4.5.25) 

This provides only an upper bound on g^^^.j{C)i rather than a bound on its absolute 
value. However it is obvious from the definition ()4.5.4p of g that g^^M^O ~ ~9l%(0^ 
so the needed lower bound follows automatically. □ 
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With Equation ()4.5.5p . Lemma 14.5.11 means that 

|Gi(C)|= 0(L'^-2) (4.5.26) 

as desired. In d < 4, this means that for sufficiently large L we have IGlI < AL'^/'^ , and 
we use the same moment generating function argument as in the previous section, we 
see that J-{he) is differentiable at /i = for all e. 
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Chapter 5 
Conclusion 

The previous sections have concluded the proof of the rounding effect for quantum 
lattice systems; that is, that first order phase transitions (and therefore, in light of 
Corollary 12.2.41 and Theorem 12.2.51 long range order) are impossible in the presence 
of direct randomness in low dimensions. This has been done by establishing a unified 
analysis of free energy fluctuations applicable to both classical and quantum systems. 

At the same time, much remains to be said about the character of the "rounded" 
phase transitions, and of the exceptional cases which have appeared in the course of 
this work. No simple statement is likely to encapsulate the situation in this context; 
certainly none can be advanced at this time. Knowledge of this area continues to grow, 
and some techniques which may be used to shed further light on it are discussed in 
Appendix [Xl 
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Appendix A 

Methods for numerical studies of random field spin 

systems 

A.l The mciximum flow representation of the Ising model ground 
state 

The RFIM at zero temperature has the considerable vh'tue that for particular finite 
field configurations the ground state can be computed easily and exactly thanks to a 
relationship with the maximum network fiow problem. 

A maximum flow problem is the following. We are given an undirected graph (that 
is, a finite collection of vertices (points), some pairs of which are connected by edges), 
with two special vertices, the source s and the sink t; each edge has a capacity, a 
finite nonnegative number which we can denote by a symmetric matrix Cij whose 
indices label the vertices. A flow is an antisymmetric matrix Fij which does not exceed 
the capacities {\Fij\ < Cij\/i^j) and which is conserved Fij = 0) except at the 
source and the sink. A maximum fiow is one which maximizes the total current from 
the source to the sink, which is given by Fgj = Fjt- This problem has been 
extensively studied by computer scientists, and there are a number of well-studied and 
efficient algorithms for solving it. Most standard implementations (for example the 
Boost Graph Library [68j) assume that the capacities are integers, which we shall see 
is inconvenient for our purposes, but this can be circumvented by rescaling and using 
very large integers. 

The relationship to the Ising model is through the related minimum cut problem: 
given the same objects as in the maximum flow problem, a cut is a choice of a division 
of the vertices into two components, one (call it S) containing the source and the other 
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(T) the sink. Each cut has a cost, which is the total of the capacities of ah edges which 
connect S to T, Ylies X^jgt ^ minimum cut is a cut which minimizes this cost 

function. Given a ferromagnetic Ising model with Hamiltonian 

we make a graph whose vertices are the sites plus a source and sink, with edges of 
capacity Jij connecting each interacting pair of sites, an edge of capacity hi connecting 
each i with positive field to the source and one with capacity —hi connecting each site 
with negative field to the sink. Then a configuration corresponds naturally to a cut 
with S being the sites with spin +1 and the source; and T-L is precisely the cost of this 
cut, so a minimum cut corresponds to a ground state [69]. 

The "max cut - min flow theorem" [70^ [7T] provides a connection between these 
two problems. It states that in a maximum flow, the saturated edges (those with 
\Fij\ = Cij) divide the graph in such a way as to provide a minimum cut (or several, if 
they divide the graph into more than two connected components), and that all minimum 
cuts for a given problem can be obtained in this way. The basic idea (also used in the 
popular push-relabel algorithm to solve the problem [72l [73]) is that if there is a path 
of unsaturated edges connecting the source to the sink then it is possible to increase 
the total current by increasing the flow along each edge of that path. 

Together with modern algorithms for solving the max flow problem, this allows the 
zero temperature random field Ising model to be simulated very efficiently [32 \ 1331 [35]. 
avoiding the extremely slow convergence which plagues monte carlo studies of disordered 
systems at low temperature. This has been the main method used for numerical studies 
of the random field Ising model, although in the last few years histogram reweighting 
methods like the Wang-Landau algorithm [341 [35] have made finite temperature simu- 
lations practical as well. 

The maximum flow representation also allows the following more general derivation 
of the exact solution of the one dimensional random field Ising model described in Sec- 
tion [TTLT] above. Construct two different graphs (as shown in Figure [AT]) by beginning 
with the construction in the previous section, only for all sites "before" x; then add an 
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Figure A.l: Example of the maximum flow graph used to derive Equation (ll.l.Sp 

edge with capacity J from site x — 1 to the sink to make the first graph, and instead 
connecting to the source to make the second graph. If the maximum flow for the first 
graph has nonzero flow through the new edge, take this as Ux] otherwise, is minus 
the flow in the new edge of the second graph. It is easy to see that this gives a unique 
set of values which satisfy Equation (|1.1.5p , and doing the same with the other half of 
the system does the same for Vx- 

We can now obtain the value of cr^ as follows. Suppose > 0; then we will have 
<Tx = 1 if there is a maximum flow in which the corresponding edge does not saturate 
{Fsx < hx). We can obtain a maximum flow for the whole system by pasting together 
the flow graphs representing the different parts as used above to obtain Ux and Vx', then 
graph can accommodate a flow Fgx of up to —Ux — Vx, but no larger; so hx > —Ux — Vx 
implies ax = 1, and hx < —Ux — Vx implies ax = —1. The same follows by a similar 
examination of the cases hx < and hx = 0. 

A. 2 Monte Carlo methods for XY and clock models 

A. 2.1 The lookup table algorithm for the Clock model 

As noted in Section ll.l.2| the clock model has been frequently used as a substitute for 
the XY model for reasons of computational efficiency. The metropolis algorithm has 
a high rejection rate at low temperatures. If one tries to avoid this by using the heat 
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bath algorithm, it appears that one has the choice of samphng the distribution of trial 
moves with either rejection sampling (which in effect reproduces the same problem) or 
through an inverse transform method (which is slowed down by the need to evaluate 
a large number of transcendental functions). Heat bath updating in the clock model 
can be implemented much more rapidly by compiling a lookup table, since only a finite 
number transition probabilities need to be calculated for a given set of parameter values. 
A new table must be calculated whenever one changes the random field distribution, 
coupling constant, or temperature, but can be reused for different configurations of the 
random field. 

To be more precise, if we rewrite the Hamiltonian (jl.l.Sp as 

/3n = -Y,ij Yl ^y + hA-a,, (A.2.1) 

X \ \y-x\=l / 

then range of relative energies involved in rotating a single spin are controlled by an 
effective field kx ■= J '^\y-x\=i ^y + ^x- If we denote the number of allowed spin values 
by g, and use a random field taking nq values with the same symmetry, then thanks to 
the symmetry in interchanging the neighboring spins it takes no more than nq^'^ values, 
which are related by a g-fold symmetry; for each of these a table of g — 1 elements needs 
to be recorded to specify the transition probabilities, giving a table with 

T = nq^'^-^{q-l) (A.2.2) 

elements, usually 32-bit integers (there is a redundancy in this description, and in fact 
kx takes no more than nq x (^'^^^~^) values, but taking full advantage of this complicates 
the algorithm). For g = 12, d = 3, n = 2 (the most ambitious case I know to have been 
implemented [l3]), this gives a table with about 6.6 x 10^ elements. Once this table 
has been calculated, each update step involves only a small number of arithmetical or 
logical calculations and the generation of a single pseudorandom number. Assuming 
that the table is stored in random access memory so that the time required to retrieve 
a specified element is independent of the table size, we can examine the computing 
time required by dividing the algorithm used to sample the update distribution into 
the following steps: 
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1. Determine which value of k to use; at most d+1 steps of constant complexity 

2. Generate a pseudorandom integer r 

3. Successively look up the probability p of the candidate configurations; if p > r, 
choose that configuration; otherwise set r ^ r — p and move on. This is at most 
q steps of constant complexity. 

From this, we can confidently expect that the computer time required for each monte 
carlo step should not grow faster than linearly in q so long as the table fits in available 
random access memory. 

A. 2. 2 A modified Ziggurat algorithm for the XY model 

I will first describe a slightly modified Ziggurat algorithm for a random variable taking 
values on [0,7r] with a decreasing probability density function p{x), before moving on 
to a further modification to accommodate the situation relevant to the XY model. In 
a preparation step, one approximates the graph of p{x) with a collection of N boxes 
(see Figure of height hi width Wi, and left coordinate x^; there is an easy method 
for choosing these parameters so that the boxes have equal volume. To include some 
boundary cases, we take xn+i = vr and /iat+i = p(vr). 

One can then sample the desired distribution in the following way: 

1. Choose one of the boxes i at random with equal probability 

2. Independently and uniformly choose a random number x from [xjjXj-i-i] and y 
from [0, hi]. If y < /li+i, return x. 

3. Otherwise, calculate p{x). If y < p{x) return x, otherwise start over with step 1. 

This is a clever way of doing rejection sampling: uniformly select a point in the union of 
the boxes, accept it if it is under the graph of the desired probability density, otherwise 
reject it. When is reasonably large, to begin with the points generated will be 
accepted most of the time, and in addition they can usually be accepted without even 
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Figure A. 2: The ziggurat algorithm: a probability density, bounding boxes of equal 
volume 

computing p{x), which can lead to sampling which is even more efficient than an inverse 
transform method. 

In heat bath simulations of the XY model, the key (unnormalized) probability 
distribution is 

Pk{0) =e^p{k cos 9), (A. 2. 3) 

where k is the magnitude of the local field described above, and 9 is the angle between 
that field and the new spin direction. This is monotone on [0, vr] and one can extend 
it to the full range by adding a step which reflects the spin with probability 0.5. The 
problem with using the Ziggurat method is that the probability distribution depends 
on a continuous parameter k, however it can be generalized in the following manner to 
accommodate this situation. 

The idea is illustrated in Figure lA.31 For a given range [ki,k2], generate a set of 
boxes as in the original Ziggurat algorithm, but with maxfcj<fc<fc2 Pa:(2;). In addition, 
for each box we also need to record ti = mmkj^<ck<k2 Pki^i+i) ■ Then we can sample pk 
as follows: 
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Figure A. 3: The modified ziggurat algorithm: two probabihty densities, bounding boxes 
of equal volume. The probability densities shown are those of Equation ()A.2.3p . with 
k = I and k = 1.5, rescaled for pfc(O) = 1. 

1. Choose one of the boxes i at random with equal probability 

2. Independently and uniformly choose a random number x from [a;j,Xj-|-i] and y 
from [0, hi]. If y < tj, return x. 

3. Otherwise, calculate Pk{x). If y < Pk{x) return x, otherwise start over with step 1. 

The efficiency of this method depends on whether we can partition the relevant range 
of k so that the variation of pk is small enough that the rejection rate and the frequency 
with which p^ is calculated do not increase to much. If we denote the maximum value 
of the random field by H, then k runs from zero to 2Jd + H. p^ is monotone in k 
(although whether it is increasing or decreasing depends on x and the normalization 
used) which makes calculating minima and maxima with respect to k very simple. For 
the ranges of k relevant to simulating 3-dimensional systems near the apparent critical 
temperature (roughly J = 2), it is feasible on a computer with 2GB of RAM to store a 
table which requires calculating p^ less than one time in 1000. 
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Appendix B 

Rounding of First Order Transitions in Low-Dimensional 
Quantum Systems with Quenched Disorder 

(With M. Aizenman and J. L. Lebowitz. Published as [1].) 
Abstract 

We prove that the addition of an arbitrarily small random perturbation to a quantum 
spin system rounds a first order phase transition in the conjugate order parameter in 
d < 2 dimensions, or for cases involving the breaking of a continuous symmetry in 
d < 4. This establishes rigorously for quantum systems the existence of the Imry-Ma 
phenomenon which for classical systems was proven by Aizenman and Wehr. 

A first order phase transition, in Ehrenfest's terminology, is one associated with a 
discontinuity in the density of an extensive quantity. In thermodynamic terms this 
corresponds to a discontinuity in the derivative of the free energy with respect to one of 
the parameters in the Hamiltonian, more specifically the one conjugate to the order pa- 
rameter, e.g. the magnetic field in a ferromagnetic spin system. In what is known as the 
Imry-Ma phenomenon [H [11] , any such discontinuity is rounded off in low dimensions 
when the Hamiltonian of a homogeneous system is modified through the incorporation 
of an arbitrarily weak random term, corresponding to quenched local disorder, in the 
field conjugate to the order parameter. 

This phenomenon has been rigorously established for classical systems j251[3]. where 
it occurs in dimensions d < 2, and d < 4 when the discontinuity is associated with 
the breaking of a continuous symmetry. In this letter we prove analogous results for 
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quantum systems at both positive and zero temperatures (ground states). 

The existence of this effect was first argued for random fields by Imry and Ma on the 
basis of a heuristic analysis of free energy fluctuations. While the sufficiency of Imry 
and Ma's reasoning was called into question, the predicted phenomenon was established 
rigorously through a number of works [ini Ell 1221 ESI EH E51 E]- The statement was 
further extended to different disorder types by Hui and Berker [26\ [71| . 

The general existence of the Imry-Ma phenomenon in quantum systems was not 
addressed by these rigorous analysis, and in particular the Aizenman-Wehr [251 [3] proof 
of the rounding effect applies only for classical systems. However, as stressed in [55] . 
establishing whether the Imry-Ma phenomenon extends to first order quantum phase 
transitions (QPTi) is an important open problem. The results presented here answer 
this question. We find that the critical dimensions for the phenomenon for quantum 
systems are the same as for classical systems, including at zero temperature. 

We consider spin systems on the d-dimensional lattice Z*^, where the configuration 
at each site is described by a finite-dimensional Hilbert space, with a Hamiltonian of 
the form 



where {k^} are translates of some local operator kq, and h and e are real parameters. 
The quenched disorder is represented by {ijx}, a family of independent, identically 
distributed random variables. T-Lq may be translation invariant and nonrandom, or it 
can include additional random terms (although we will not discuss the latter case, our 
results hold there also). For convenience we will assume that \\kx\\ = 1, which can be 
arranged by rescaling h and e. We will refer to the r]s as random fields, although in 
general they may also be associated with some other parameters, e.g. random bond 
strengths. 

An example of a system of this type (with Kx = (Tx ) is the ferromagnetic transverse- 
field Ising model with a random longitudinal field [53] (henceforth QRFIM), with 



where {i = 1,2,3) are single-site Pauli matrices, and Jx-y > 0. The QRFIM has 




(B.0.1) 



X 




(B.0.2) 



63 



recently been studied as a model for the behavior of LiHoxYi_xF4 with x > 0.5 in a 
strong transverse magnetic field \75\ [76] . 

We will examine phase transitions where the order parameter is the volume average 
of the expectation value of with respect to an equilibrium (KMS) state, and show 
that this quantity cannot be discontinuous in h for low-dimensional systems. As is well 
known, this order parameter is related to the directional derivatives (it) of the free 
energy density, 

m±{T,h,e) := -^HT,h,e) (B.0.3) 
where, as usual, at positive temperatures 

T(T,h,e)= lim -^logTre-'^^r (B.0.4) 

(with /3 := l/ksT), and T{0, h, e) is the corresponding limit of the ground state energy. 
Here Hr is the Hamiltonian of the system restricted to the finite box F C Z"^, and |r| is 
the number of sites in that box. It is known under the assumptions enumerated below 
that for almost all rj this limit exists and is given by a non-random function of the 
parameters (see, e.g. [3l[62]), which does not depend on the boundary conditions. By 
general arguments which are valid for both classical and quantum systems, J- is convex 
in h; therefore the directional derivatives exist, and are equal for all but countably 
many values of h }46j . 

For typical realizations of the random field, the interval [m^{T,h,e),m^{T,h,e)] 
provides the asymptotic range of values of the order parameter for any sequence of 
finite volume Gibbs states or ground states (the argument is similar to that found in [3] 
for classical systems). At a first order phase transition m_ < m+, and there are then 
at least two distinct infinite volume KMS states |46| with different values of the order 
parameter. In the QRFIM the m+ and the m_states can be obtained through the -|- or 

(3) 

— boundary conditions (i.e. the spins a-h are replaced by ±1 for all x ^ F). In general, 
such states are obtained by adding ±5 to the uniform field h and letting (5—7-0 after 
taking the infinite volume limit. 

Our discussion is restricted to systems satisfying: 
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A. The interactions are short range, in the sense that for any finite box A G Z'^ the 
Hamiltonian may be decomposed as: H = H/^ + Va + H\c, with H\ acting only 
in A, H\c only in the complement A"^, and V\ of norm bounded by the size of the 
boundary: 

II^aII < C\dA\ . (B.0.5) 

B. The variables ijx have an absolutely continuous distribution with respect to the 
Lebesgue measure (i.e. one with a probability density with no delta functions), and 
a finite rth moment, for some r > 2. 

Our main results are summarized in the following two statements. The first applies 
regardless of whether the order parameter is related to any symmetry breaking. 

Theorem B.0.1. In dimensions d < 2, any system of the form of (jB.O.ip satisfying 
the above assumptions has m^{T,h,e) = m-{T,h,e) for all h, and T > 0, provided 

The next result is formulated for situations where the the first order phase transition 
would represent continuous symmetry breaking. An example is the 0{N) model with 

"Ho = - ^ Jx-y^x ■ (B.0.6) 

where a are the usual quantum spin operators. More generally, "Hq is assumed to be 
a sum of finite range terms which are invariant under the global action of the rotation 
group SO{N), and (fx is a collection of operators of norm one which transform as the 
components of a vector under rotations. With the random terms the Hamiltonian is 

n = 'Ho-Y,ih + effx)-ax. (B.0.7) 

Theorem B.0.2. For the SO{N)- symmetric system described above, with the random 
fields fjx having a rotation-invariant distribution, the free energy is continuously differ- 
entiable in h at h = whenever e / 0, d < 4, and N > 2. 

Before describing the proof, let us comment on the implications of the statements, 
and their limitations. 
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1. While the statements estabhsh uniqueness of the expectation value of the bulk 
averages of the observables Kx, or ffx (in Theorem IB.0.2p . they do not rule out the 
possibility of the coexistence of a number of equilibrium states, which differ from each 
other in some other way than the mean density of k, which they share. More can be 
said for models for which it is known by other means that non-uniquess of state is 
possible only if there is long range order in k. (Such is the case for QRFIM, through 
its relation to the classical ferromagnetic Ising model in d + 1 dimensions [52j.) 

2. The results address only the discontinuity, or symmetry breaking (as in the QR- 
FIM), but they leave room for other phase transitions, or singular dependence on h. 
For instance, for the Ashkin- Teller spin chain for which Goswami et. al. [55] report 
finding the Imry-Ma phenomenon in some range of the parameters but not elsewhere, 
the results presented here rule out the persistence of a first-order transition between the 
paramagnetic and Baxter phases in the full range in the model's parameters. However, 
they do not rule out the possibility of other phase transitions. 

3. Randomness which does not couple to the order parameter of the transition need 
not cause a rounding effect. For example, in the transverse-field Ising model in a random 
transverse field, where the random field r]x couples to (t^^\ ferromagnetic ordering 
is known to persist [53^ I52j . Presumably the same is true for the Baxter phase of 
the Ashkin- Teller model. It was even suggested that there are systems in which the 
introduction of randomness of this sort may even induce long range order which would 
not otherwise be present [771 [78], ^-iid our results do not contradict this. In addition 
we can draw no conclusions about quasi-long-range order, that is power law decay of 
correlations, including of Kx- 

Other comments, on the technical assumptions under which the statements hold, 
are found after the proofs. 

The proofs of Theorems IB . . II and IB . ."2] are based on the analysis of the differences, 
between the m-|_ and the m_-states, in the free energy (at T = 0, ground state energy) 
which can be ascribed to the random field within a finite region A of diameter L. 
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Putting momentarily aside the question of existence of hmits, a relevant quantity could 
be provided by: 



lim lim Av 



Av 



A.r 



(B.0.8) 



where p is the difference of free energies 



with 



-1 



T •- ^logTr exp(-/3i/^''') , 



(B.0.9) 



(B.0.10) 



T]^^^ is the random field configuration obtained from rj by setting it to zero within A, 
and Av[-|r/A] is a conditional expectation, i.e. an average over the fields outside of A. 
(The modification of the field h by ±(5 serves to select the desired {m±) states). 

Somewhat inconveniently, it is not obvious that for all models the limits in (lB.0.9h 
exist. Nevertheless, one can prove that for each system of the class considered here there 
is a sequence of volumes Tj "L^ for which the limit exists for all A, with convergence 
uniform in r^A- The proof of this assertion is by a compactness argument, whose details 
can be found elsewhere [79 1. 

The essence of the proof of Theorem IB. 0.11 is the contradiction between two esti- 
mates: 



i. Under Assumption lAl equation (jB.O.SP : 



|GA(r/)| < 4C|M| 



(B.0.11) 



ii. Whenever m_ < ?ti+, Ga/\/|A| converges in distribution to a normal random 
variable with a positive variance (as one would guess by considering the difference in 
the random field terms between states of different mean magnetizations, neglecting 
the states' local adjustments to the random fields). 

More explicitly, for the upper bound we note that in the absence of the interaction 
terms Va, the right hand side of (|B.0.9p would be zero. Using (|B.0.5p . one gets (|B.0.1ip . 



^ Sections lO and lO 
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To prove the normal distribution for G^p, we apply a theorem of [3] (Proposi- 
tion 6.1) (as corrected in [56], p. 124). It implies that for A Z"^, under Assumption [Bl 
Ga/y^IAI converges in distribution to a normal random variable with variance of the 
order of 



Av 



7Ti_|_ — m_ 



(B.0.12) 



The two statements described above contradict the assumption that m_ < in 
dimensions d < 2. That is so even at the critical dimension, where L*^/^ = L'^~^. The 
reason is that the lower bound implies the existence of arbitrarily large fluctuations on 
that scale, whereas the upper bound is with a uniform constant. This proves Theo- 
rem EHU 



The above proof is similar to that of the classical results [3l [25] which this work 
extends. However, the discussion of the free energy fluctuations was based there on the 
analysis of the Gibbs states, and more specifically of the response to the fluctuating 
fields of the 'metastates' which were specially constructed for that purpose. Except 
for special cases, such as the QRFIM, that argument was not available for quantum 
systems, where the equilibrium expectation values are no longer given by integrals over 
positive measures. The proof of the quantum case is enabled by a more direct analysis 
of the free energy. 



Theorem lB.0.2l is proven by establishing that in the presence of continuous symmetry 
the upper bound (jB.O.lip . for A = [— L, L]"', can be replaced by: 

\GA{riA)\<KL''-\ (B.0.13) 

Here G\ is defined as in (|B.0.9p . (|B.0.8p . but h = 0, and 5 is replaced by by (5 := 6e with 
e a unit vector. This change in the upper bound raises the critical dimension to d = 4. 
To obtain (|B.0.13P we focus on 



aiAvA) ■= Av 



5e 



F^ 



VA 



(B.0.14) 



Since GA{r]A) = lim^^^o hmr-+zd | (5A,r(^A) " di^ri^)) ' any uniform bound on \g{j.\ 
for given A implies a similar bound on IGaI- The claimed bound may be obtained 
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though a soft-mode deformation analysis, which we shall make explicit for the case of 
pair interaction (the general case can be treated by similar estimates). 

The free enrgy F^'^^^ in ()B.0.14p may be rewritten by rotating both the spins and 
the field vectors with respect to an axis perpendicular to e at the slowly varying angles 

||x|| < L 
-j^ — ^vr, L < \\x\\ < 2L ■ 



0, 

\\x\\-L^ 



(B.0.15) 

vr, ||x|| > 2L 

The rotation aligns the external fields in the two terms (±(5e), except within A where 
the effect is negligible when 5 — )• 0. The effect of the rotation on the random fields 
is absorbed by rotation invariance of the average. In the end, the Hamiltonian of the 
rotated system differs from the Hamiltonian used to define the other free energy by 



A-f/e := ^ Jx-y [ 

{x,y}cr 



(B.0.16) 



When the resulting expression for F^'~^^ in ()B.0.14p is expanded in powers of 9x — 
6y ~ ttWx — y\\/L, the zeroth-order term cancels with Fy~^^, and the second and higher 
order terms yield the claimed bound. The main difficulty is to eliminate the first order 
terms, which amount to a sum of 0{L'^) quantities each of order 1/L. However, the 
sign of these terms is reversed when the rotation is in the reversed direction. To take 
advantage of this, we combine two expressions for 5^p(?/a) with the rotations applied 
in opposite directions, yielding: 

5A,r(^A) = Av [logTre-^^H 

-i logTr e-/'(^+^^^) - i log TV e-^(^+^^-^) 



(B.0.17) 



(where H = Hll' By known operator inequalities [46]: 



logTr e-^^-UogTv e-/^(^+^^£) - ilogTY e-/^(^+^^-^) 



(B.0.18) 



The right hand side is zero to first order, and one is left with an upper bound on 17^ p 
of the desired form. Repeating this analysis with the roles of the terms exchanged we 
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obtain an identical lower bound, and thus inequality (jB.O.lSP follows. 

The above argument is spelled out in detail in [7^^. Let us end with few additional 
comments on the assumptions. 

4. For Theorem IB. 0.21 the assumption that the interaction has a strictly finite range 
can be weakened to a condition somewhat similar to Assumption [Xj For pair interac- 
tions (equation ()B.0.7P ) it suffices to assume: 



5. The restriction to absolutely continuous distribution excludes a number of models 
of interest. Such an assumption is generally necessary at zero temperature, as can be 
seen by the behavior of the Ising chain in a random field [28J which takes only a finite 
number of values. For positive temperatures it can be replaced by the requirement that 
the distribution has a continuous part which extends along the entire range of values. 
For the QRFIM at finite temperature one need only assume that the random field has 
more than one possible value, and this may well be the case more generally. 
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(B.0.19) 



^See also Section l475l 



70 



Appendix C 

Some concepts and results in mathematical probability 

In this appendix, I will review some of the probabilistic terminology and concepts used 
in this dissertation for the benefit of those readers who may find them obscure. Proofs 
and further details can be found in standard probability textbooks such as |66j or |80j . 

C.l (j-algebras and measures 

We begin with some set X, for example the set of outcomes of a class of measurements. 
Assigning probabilities to members of X is a problem when X is uncountably infinite 
(e.g. when it is a real interval), so we will instead assign probabilities to subsets of X 
(events). It turns out to be impossible to consistently do this for all subsets of X what- 
soever, but we at least want to be able to look at complements of sets (corresponding to 
logical not) and countable unions (corresponding to logical and), intersections (logical 
or) then comes for free. A collection of X of subsets of X which includes X and is 
closed under complement and countable union is called a a-algebra. The most common 
example is the Borel algebra on the real line, which is the smallest cr-algebra which 
includes all intervals. The combination of a set and a u-algebra on that set is called a 
measurable space, and in this context an element of the u-algebra is called a measurable 
set. 

A function /u assigning a nonnegative real number (possibly oo) to each member of 
a (7-algebra is called a measure if /i(0) = and 



\i=i / i=i 

for any countable collection of disjoint (nonoverlapping) sets Ai ^ X . /i is called a 
probability measure if ii{X) = 1. This definition immediately implies a number of 




(C.1.1) 



71 



the usual rules of probability: probability of mutually exclusive events is additive, the 
probabilities of a complete set of mutual exclusive events is 1, etc. A measurable space 
with an associated measure is a measure space, or if that measure is a probability 
measure it is a probability space. 

One measure we use frequently is the Lebesgue measure A, also known as the uniform 
measure on the real line. This has the property that for any finite interval A(/) is the 
length of /, and it is the only measure on the Borel algebra with this property. Another 
important Borel measure is the Dirac measure, defined by 

1, X e A 



(C.1.2) 

0, x^A 



It is possible to define integration with respect to a measure, but only with respect to 
measurable functions. A function / from one measurable space to another is measurable 
if the preimage of a measurable set is a measurable set (that is, if the set of x which 
produces some measurable range of outcomes f{x) is guaranteed to be measurable). 
This means that knowing which measurable sets contain x is enough information to 
deduce the value of /(x). Continuous functions are always measurable with respect to 
the Borel algebra. Given a measure, it is possible to define integration of a function / 
with respect to a measure mu, denoted 

J fd/j (C.1.3) 

first for positive measurable functions, then for functions whose positive and negative 
parts give a finite integral (these are then integrable functions). This has many of the 
usual properties of an integral. We can also define integrals over a measurable set A, 
which have the property that 

dfi = /i(y4). (C.1.4) 



I A 

We can also write this in terms of the indicator function 



as 



I[A]{x):=\ ^' (C.1.5) 
[ 0, x^A 



Jjdf^ = J f^i^^df, (C.1.6) 
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Integration with respect to the Lebesgue measure (leaving aside some issues of 
divergence) coincides with the usual notion of integration on the real line: 

/oo 

f{x)dx (C.1.7) 
-oo 

fdX= j f{x)dx (C.1.8) 

[a,b] J a 

whenever the relevant expressions are all well defined. A set A which is contained in a 
measurable set B with /i(-B) = is called a null seqj; two functions which differ only 
on a null set have the same integrals. 

Given a measure jjL and a nonnegative measurable function /, we can define a mea- 
sure by 

f^x{A) = [ fdpi. (C.1.9) 

J A 

If it is possible to write a measure v as //u, then we saw that v is absolutely continuous 
with respect to /i; / is the Radon-Nikodym derivative or density of with respect to 
fi. When I say that a measure is absolutely continuous without specifying another 
measure, this should be understood to be the Lebesgue measure. 

C.2 L'P norms and spaces; convergence of measurable functions 

Consider some measure space {X, X, fi), and the following functional on the measurable 
functions (1 < p < oo): 

i/r^^j (c.2.1) 

This is not quite a norm, since it is unaffected by changing the value of / on a null 
set. It is however a norm on the resulting equivalence classes, and so it is frequently 
referred to as the norm; the space of equivalence classes of functions with ||/||p finite 
is denoted C^ip.)- This can be extended to p = c« by setting 

oo = , if sup|/(x)|. (C.2.2) 

A:fi{X\A)=0 x&A 



^This is not precisely the same as saying that A has measure zero, since A is not assumed to be 
measurable; however this distinction is rarely important. 
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For any 1 < p < oo, the norm provides a notion of convergence like any other norm. 
We say that a sequence /„ converges to / in >C^(/u) iff 

lim =0. (C.2.3) 

n— >-oo ^ 

One useful property of these norms is H aider's inequality, which states that 

Il/fflli<ll/llpll5ll, (C.2.4) 

whenever l/p + l/q = 1 (including p =l,q = oo). 

There are several other important notions of convergence for measurable functions. 
The first is convergence in measure: fn^fi^^ measure if 

l^{{\f-fn\>e}nA)^0 (C.2.5) 

for any e > and any measurable set A with fi{A) < oo. 

A sequence /„ converges to / almost everywhere (in a probability space, almost 
certainly or almost surely) if there is a null set N so that 

fn{x) -> /(x)Vx ^ N. (C.2.6) 

Almost everywhere convergence and convergence both imply convergence in measure. 

The following results give some relationships between convergence of functions and 
of integrals. 

Theorem C.2.1 (Beppo-Levi) . Let fn be a sequence of integrable functions with 

j fndfj, < OO, with fn f (monotone convergence) almost everywhere, with f mea- 
surable. Then 

lim / fnd^i = f fdn. (C.2.7) 

Theorem C.2.2 (Lebesgue, dominated convergence). Let fn be a sequence in L^{ii) 
with fn ^ f i'f' measure for some measurable f . Assume there is a nonnegative 
g G L^{fi) such that |/n| < 5 almost everywhere for all n. 
Then f G L^{ii) and fn^ f in '> in particular 

j fndfi ^ j fdfi. (C.2.8) 
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We can also talk about limits of sequences of measures on a metric space (e.g. R). 
The most common notion is the following: a sequence of finite measures /x„ on the same 

measurable space converges weakly to ii if 

J fdfin ^ J fdfi (C.2.9) 
for all bounded, continuous functions /. 

C.3 Random variables, expectations, conditional expectations 

An integrable function on a probability space is called a random variable. The average 
(or expectation) of a random variable (j) is 




(C.3.1) 



The expectation satisfies Chebyshev 's inequality, 

Av I[\(p -Av(P\>e]< e-^ {Av{(j)'^) - (Av ^f) . (C.3.2) 

It is possible to say a great deal about a random variable without specifying which 
probability space it is defined on. The important fact is that a random variable (p 
defines a probability measure, called its distribution and denoted by P,^, on the space 
in which it takes its values. This is exactly the probability that 4> will take a value 
in a particular measurable set. We can talk about multiple random variables being 
identically distributed if their distributions are the same. Random variables converge 
in distribution if their distributions converge weakly. One important distribution is the 
normal distribution with mean m and variance b^, denoted N{m, 6^), which is absolutely 
continuous with respect to the Lebesgue measure and has density e~^^/^''^/(6\/27r). 

Among the important properties of a random variable are its moments. The kth 
moment is Av/^'. A bounded random variable is characterized by its (integer) moments, 
or by its exponential moment generating functional Av e^^ . 

Random variables (j) and ip defined on the same probability space are independent 
or independently distributed if, for any measurable A and B, the probability that (f) E A 
and '0 G -B is the product of the two separate probabilities. A similar definition holds for 
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arbitrary collections of random variables. We frequently speak of collections of random 
variables being independently and identically distributed, abbreviated i.i.d.. 

Given a random variable (j) and a sub-ci-algebra of the one for which it was defined, 
the conditional expectation of with respect to (written Av [i?!*]-?^]) is another random 
variable which is measurable with respect to J- and which satisfies 

Av {(t>I[A]) = Av (Av [01^] I[A]) (C.3.3) 

for any A £ J^. A version of Theorem IC. 2. 21 holds for conditional expectations: 

Theorem C.3.1. Let phin be a sequence of random variables, and let be a nonneg- 
ative random variable on the same space, with {(pnl < o,nd cpn ^ P almost surely. 
Then 

lim Av [(l)n\T] = Av (C.3.4) 

n— >oo 

almost surely and in L^. 

There is a great deal to be said about collections of random variables which fail to 
be independent in particular ways. To talk about this, we first introduce the idea of a 
filtration, which is a nested sequence of c-algebras, J-iCJ-2C---. A stochastic process 
is a sequence of random variables such that each is measurable with respect to 
the corresponding element of the filtration. 

A martingale is a stochastic process which satisfies 

Av [(pnl^m] = Av [(pmlJ^m] (Vn > m). (C.3.5) 

There are a number of results about the limits of martingales, which are referred to at 
various points in the text. 
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Appendix D 

Product Measure Steady States of Generalized Zero 

Range Processes 

(With J. L. Lebowitz, published as |5]) 
Abstract 

We estabhsh necessary and sufficient conditions for the existence of factorizable steady 
states of the Generahzed Zero Range Process on a periodic or infinite lattice. This 
process allows transitions from a site i to a site i + q involving (a bounded number 
of) multiple particles with rates depending on the content of the site i, the direction 
q of movement, and the number of particles moving. We also show the sufficiency of 
a similar condition for the continuous time Mass Transport Process, where the mass 
at each site and the amount transferred in each transition are continuous variables; we 
conjecture that this is also a necessary condition. 

D.l Introduction 

The classical zero range process (ZRP) is a widely studied lattice model with stochastic 
time evolution [81]. To define the process consider a cubic box A C with periodic 
boundary conditions, i.e. a d-dimensional torus. At each site z of A there is a random 
integer-valued variable iii £ {0,1,...}, representing the number of particles at site i. 
The time evolution is specified by a function aq{ni) giving the rate at which a particle 
from a site i containing rii particles jumps to the site i + q, where q runs over a set of 
neighbors E (the most common choice is E = {=bei, ±62, . . .}, but our treatment holds 
for any finite E which spans Z*^). The name zero range indicates the fact that the jump 
rate from i to i + q depends only on the number of particles at i. 
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It is easy to see, when the system is finite, aq{n) + a_q(n) > S > for all re > 
(we always have a(0) = 0) and E spans Z,'^, that all configurations with a given total 
particle number = X^ieA ''^i mutually accessible, and hence there is a unique 
stationary measure -Pa(ii;-^) for each N. Normalized superpositions of these measures 
yield all of the stationary states of this system. Conversely, given a stationary measure 
for which there is a nonzero probability of N particles being present in the system one 
can obtain Pa{e.',]^) by restricting that measure to configurations with N particles. 

The ZRP was first introduced in [82j . It was assumed there and in most subsequent 
works that the rates ag(re) are of the form 

ag{n) = gqa{n) (D.1.1) 

with Qq independent of n and a(n) independent of q. In this case the system has the 
unique steady state given by [82l [83] 

VieA / ieA 



where 



P^"^^ = Tvi 7u\ (D.1.3) 

Cm and c are normalization constants given by 

The unique stationary measure Pa (re; N) is thus a restriction to configurations with N 
particles of the product measure rijeA^'('^0 with single-site distribution p{n). 

In the limit A — )■ Z"' with A^/|A| — t- p the only stationary extremal measures, i.e. 
the only stationary measures with a decay of correlations, are product measures with 
pin) given in Equation (ID.l.Sp as the distribution of single site ocupation numbers 
|83j . These states are parameterized by A, which plays the role of the fugacity in 
an equilibrium system, with different values of A corresponding to different expected 
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particle densities p, where 



oo oo 



nX 



(D.1.6) 



n=l n=l 



Recently there has been a revival of interest in the ZRP. For certain choices of a(n), 
for example when a{n) ~ 1 + h/n for large n, the ZRP on Z exhibits a transition 
between a phase where all sites almost certainly contain finite numbers of particles to a 
'condensed' phase where there is a single site containing an infinite number of particles 
|84^ 185], This condensation has attracted attention as a representative of an interesting 
class of phase transitions in one dimensional non-equilibrium systems, and has also been 
applied to models of growing networks [811 [86] . 

Evans, Majumdar and Zia [87] have proposed a generalization of the ZRP, called 
a Mass Transport Model (MTM). They considered a one dimensional lattice on which 
there is a continuous 'mass' m-j > at each site, with a parallel update scheme in which 
at each time step a random mass /ij, < /Uj < m-j moves from each site i E Z to the 
neighboring site i + 1 with a probability density (j){fj,,mi). This process shares many 
features with the (totally asymmetric) Zero Range Process, in particular the existence 
of a condensation transition in certain cases [88]. One very significant difference from 
the ZRP, however, is that the system has a product measure steady state if and only 
if m) satisfies a certain condition. Taking a limit in which the probability of a 
transition at any given site at any given time step goes to zero (discussed in [87|) gives 
a stochastic process on continuous time (equivalent to a model with random sequential 
local updates) which we will call a Mass Transport Process or MTP, in which the rate of 
transitions from site i to i + 1 is given by a(/i, m). This process has a product measure 
steady state if and only if there exist functions g and p such that 



An interesting question, then, is whether similar criteria for the existence of a prod- 
uct measure exist for such processes in higher dimension and with movement in both 
directions allowed. This question is already relevant for the ZRP. A particular case in 




(D.1.7) 
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d = 



2, studied in [891 [90], has 



a±i{n) = a[l - 6n,o] 
a±2{n) = a^'^\n) = an 



(D.1.8) 



(D.1.9) 



i.e. a constant rate (independent of n) per occupied site for moving in the ±2; direction 
and a rate proportional to n in the iby direction. A treatment of this system based 
on fluctuating hydrodynamics and computer simulations (originally conducted on a 
similar but not quite equivalent system, but which we have reproduced on this system) 
suggests that this particular system has correlations between occupation numbers at 
different sites a distance D apart decaying according to a dipole power law D-^. This 
behavior, which is very different from a product measure steady state or its projection 
()D.1.2p . is conjectured to be generic for nonequilibrium stationary states of systems 
with non-equilibrium particle conserving dynamics in d > 2. 

In the present work we prove rigorously that Equation (jP.l.ip is a necessary and 
sufficient condition for the existence of product measure steady states for ZRPs, and 
as a consequence that the system described by ()D.1.9p has no product measure steady 
states. This condition in turn is a special case of a condition on a class of systems which 
we call Generalized Zero Range Processes (GZRP), in which we also allow transitions 
in which more than one particle moves at a time, although we will assume that the 
number of particles moving in a single transition is bounded. The rate now depends 
on the number of particles i' which move in the transition as well as the number of 
particles n at the site before the transition, and so the rates are given by a function 
aq{iy, n) with some i^max such that aq{n, z^) = whenever 1/ > Vmax- The classical ZRPs 
discussed above are a special case with Vmax = 1 and ag(l, n) = aq{n). We prove that a 
necessary and sufficient condition for the GZRP to have product measure steady states 
is 



for some non-negative gq{v) and f{n) with ^ f{n) < 00. This has a clear similarity to 
(ID.1.7D . and when i^rnax = 1 (|D. 1.1011 reduces to (ID.l.lD . 

We will prove Equation (|D.1.10p in the course of finding a weaker result for the 




(D.1.10) 
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continuous-time Mass Transport Process generalized to dimension d > 1 and to transi- 
tions in all directions. We show that these systems have product measure steady states 
when aq{fi,m) = gq{fi)p{m — fi)/p{m). This condition is also necessary under certain 
conditions (generalizing ()D.1.7p to higher dimension) and we conjecture that this is so 
in all cases. 

D.2 Factorizability in the Mass Transport Process 

Let Pkink^t) be the time-dependent probability density of finding the system in a par- 
ticular configuration m with mass at site i G A, rrii £ (0, oo). As noted above, we are 
considering periodic boundary conditions; this case is somewhat simpler than those of 
other boundary conditions. The master equation describing the evolution of Pa (m, t) 
is 



dPA{m,t) 

at 



(f-TUi 
/ diJ,ag{fi,mi)P{rn,t) 



(D.2.1) 



where 



i.q.u 

m- 



(D.2.2) 



i^j - /i, j = i 
rrij + j = i + q 

A stationary state of the system is a distribution Pa (m) such that dP\ (m, t)/dt = 
whenever Pa {m, t) = Pa {m) , or equivalently 

rnii _ i-rrii 

VV / dtiag{fi,mi)PAim) = Y\Y\ (i/ia_<j(/x, m^+g //)Pa(?21*'^''') (D.2.3) 

iSA q ^eA q 

We wish to find conditions under which there is a Pa which is factorizable, that is 
which takes the form of 

PA{m) = Ylp{mi) (D.2 A) 

«6A 

Assuming that aq{^, m) > for all m > and < ^ < m for at least one q of each 
pair of opposite directions, the system in a finite torus A will have a unique steady state 
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corresponding to each value of the total mass M = 'Y^rrii. Any linear combination of 
such states is also a solution of ()D.2.3p . Given a factorizable steady state, states of 
definite total mass can be obtained by projecting P onto the set of configurations with 
a particular value of M in analogy with Equation ()D.1.2p . 

Assume that there is a factorizable steady state as in (|D.2.4p . Let p[s) be the 
Laplace transform of p{m), and let 

(pq{fi,s) = [l/p{s)] dme~''"^aq{n,m + n)p{'m + fi) (D.2.5) 
Jo 

Note that, since aq{iJ,,m) = for m < fi, 
dme~'^"^aq{fi, m)p{m) 

= e"''^ / dnie~''"^aq{fi, m + ij.)p{m + n) = e~''^4>q{fi, s)p{s) 



We also have 

oo />m POO />oo 

dme'^"^ / diiaq{^,m)p{m) = / d^ I dme~^"^aq{fi,m)p{'m) 

^0 ^0 Jfi 

poo poo 

/ dfie'"^^ / dme'^^^aqdi, m + ix)p{m + //) 
^0 Jo 

df^e~''^'4>q{fi,s)p{s) 

(D.2.7) 

Multiplying both sides of (|D.2.3p by Y\ - e~^^"^^ and integrating over all rrii, we obtain 



oo 







TTp(sj) / dmi / dfj.aq{fi,mi)p{mi)e~''''^' 
ieA,,eE J ^0 Jo 

( \ /■°° f°° (D.2.8) 

= > 11 P(^j) / dnii / dnii+q / dfj, 

i^A^q^E J ^0 Jo 

X a-qifi, rui+q + n)p{mi - li)p{mi+q + ^)e-^.™-«.+,'».+. 

Rewriting ()D.2.8P with the aid of ()D.2.6P and ()D.2.7p and canceling common factors, 
we obtain 

poo poo 

dfi<l)qifi,Si)e-'^'' = dfi<l)^qifi,Si+g)e-'^'' (D.2.9) 

Equation (jD.2.9p will be satisfied if (though not only if) 



^qifi,s)=gqifi) (D.2.10) 
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In this case Equation ()D.2.5P gives 




/ dme ''"^ag{fj,,m + n)p{m + = gg{fi) / dme ''"'p{m) (D.2.11) 




which by uniqueness of the Laplace transform gives 



p{m — fi) 
p{m) 



(D.2.12) 



Equation (|D.2.12j) is a generahzation of the comparable formula for the unidirec- 
tional case [87j. In this case and in all other cases where, for each q £ E, either 0^ = 
or a-q = and hence either (pq = or (p^q = 0, there is in Equation (|D.2.9P only one 
term which depends on each pair mi,mi^q, and in order for the equation to be satisfied 
it must depend on only one of them. This happens only if (jD.2.10p holds for that q, 
so in these cases Equation ()D.2.12p gives the only possible rates for which there is an 
invariant product measure. 

Although in general Equation ()D.2.10p is not the only way of satisfying Equation 
()D.2.9p . solutions of this equation only correspond to realizable dynamics when p and 
ttq are non-negative and normalizable; the resulting restrictions on (j)q from Equation 
(|D.2.5p are such that it seems unlikely that there are reasonable (indeed any) rates, 
other than those in ()D.2.12p . which satisfy all of these conditions. 

Dynamics for which the system has a factorizable steady state can be found by be- 
ginning with some suitable (positive and normalizable) p{m) and then defining aq{fi, m) 
via (|D.2.12p . For example let 



where is the Heaviside step function. The possible transition rates corresponding to 
P{rn) = Y[Pc{n^) are 



where gq are arbitrary non-negative integrable functions, i.e. the rates aq{fi,m) are 
independent of m as long as p, <m. 



—cm 



9{m) 



(D.2.13) 



aq{m, n) = gq{fi)e''''9{m - /i) = gq{^i)e{m - p) 



(D.2.14) 
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D.3 Reverse processes 

In this section we will show that Equation (|D.2.12p is a necessary condition for the 
existence of factorizable steady states of any MTP whose reverse process is also an 
MTP. The relevant way in which the reverse process can fail to be an MTP is that it 
can have transition rates which depend on the mass at the target site of the transition 
as well as on the mass at the site it is leaving. 

In general, given a Markov process with transition rates K{m — )• m') and stationary 
distribution P{m), the reverse process is defined by rates K*{m — )• m') given by 

K*{m^m') = ^^M-^^^ (D.3.1) 

P{m) ^ ' 

This new process is what one obtains by running the original process backwards. Conse- 
quently the reverse process has the same stationary distribution as the original process, 
and when K is translation in variant so is 

For an MTP defined by rates ag{mi,ii), K{m — m!) is equal to aq{mi,^) for 
configurations m and m' related by moving a mass fi from site i to i + q, and to 
otherwise. The reverse process is specified by the rate function a*^{fi,m) which is the 
rate of transitions from a configuration m in which a mass /i moves from site i to site 
i + q; these are the only transitions in this process. 

When P{m) is a product measure with single-site weights p{iTi), (ID.3.ip becomes 
* / \ / , p{mi - n)p{mi+q + fi) ^Qo^ 

""-^^^^^^ = «-.(^^^^+. + ^) pi^mMm^,) ^''•'■'^ 

Rewriting, we have 

Oii,q{t^,m)—f r = a-q[^i,mi+g + fi) — (D.3.3) 

If a* defines a mass transport process, then it must be independent of all mj for j ^ i. 
In this case both sides of ()D.3.3P are equal to some function which depends only on 
and q. This which can only be true if a satisfies Equation (|D.2.12p . in which case one 
finds that 

o^lqip^m) = a-q{fi,mi) (D.3.4) 



^If for some configurations P{rn) — 0, then one defines a new process on a configuration space 
excluding these configurations (this problem does not arise in the case under consideration). 
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D.4 Factorizability in Generalized Zero Range Processes 

With mass at each site restricted to an integer particle number nj, we can reproduce 
the analysis in the previous section up to Equation ()D.2.9p . Denoting the vector of 
occupation numbers by n, and the transition rates by aq{i', n), the stationarity condition 
corresponding to Equation (ID.2.3j) is 

"i)^A(^) + + ^)^A(li'''''")) = (D.4.1) 

iGA (jS-B ^=1 

Suppose P is factorizable, 

PA{n) = llp{ni) (D.4.2) 

where p{n) is the probability of having n particles at a given site. Then define the 
generating function (discrete Laplace transform) 

oo 

p(z) = 5^zXn) (0.4.3) 

n=0 

and let 

p{z) 

Note that (j)q{i', z) >0 for all z^, z > 0. The counterpart of Equation ()D.2.9P is then 

oo oo 

Y^iM'^'Zi)= ^<0_g(z^,Zi+g) (D.4.5) 

ieA,qeEv=l ieA,qeEu=l 

We now exploit the assumption that transitions occur only for v < Vmax- Then 
choosing some j S A and q £ E and taking the kth derivative of the above expression 
with respect to Zj and zj+g gives 

u=k L'=k 

For k = Vmax^ we have 

4""""' Vmax, Zj) + </.L'™^) {ymax, Zj+g) = (D.4.7) 

For ()D.4.7p to hold for all Zj and Zj+g, both terms on the left-hand-side must be constant, 
and thus the functions 4>±q{i'max-: •) are polynomials of degree Vmax'-, being non-negative 
they must have non-negative leading terms. Equation (|D.4.7p states that pairs of leading 
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terms of these polynomials must add up to zero and so each must be zero, and therefore 
the functions (j)q{i'max, •) are polynomials of degree at most Umax — 1 for each q. 

Now setting k = Vmax — ^ we find by the same reasoning that the functions (pq^i'max, ■) 
are polynomials of degree at most I'max — 2. Proceeding in this manner we find 

^g{u,z)=gq{u) (DAS) 

as a necessary as well as a sufficient condition for ()D.4.5P to be satisfied. Referring to 
the definition of (j), this implies that 

a,(z.,n) = <7,H^^^^ (D.4.9) 
p[n) 

is a necessary and sufficient condition for the existence of a product measure. 
In the case where Vmax = 1 and ag(l,n) = ag{n), this condition becomes 




(D.4.10) 



This is what we referred to above as the classical ZRP, with the well-known stationary 
measure [HU |8l] discussed in the introduction. 

D.5 GZRPs on infinite lattices 

In order to show that we have really found all of the factorizable steady states of this 
class of systems, it remains to be established that the conditions obtained also apply 
to an infinite lattice; that is, that there are not rates for which the resulting GZRP on 
an infinite lattice has product measure steady states while the GZRPs defined on finite 
lattices have none, or do not have the same such stationary states. 

Let P{n) be a product measure with single-site distribution p(n) which is stationary 
for rate functions a on Z"', and let A be a finite box in Z"' such that there is some io & A 
such that iQ + q £ A for all q £ E. Denote by the configuration of the system inside 
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of A, and let P{n/iJ be the marginal distribution of this configuration. Then we have 
d 



rii + q 



dt 



ieA qeE iy=l i<aAqeEiy=l 
oo oo 

~E E EE^^^'^'^^^^^aX'^) 

iedAqeE:i+qeAi^=ln=u (D.5.1) 
oo n 

+ E E EE«'?(^'")^K'''M^) 

iedA q€E:i+q€A 71=1 u=l 

=0 

where dA = {i A\{3q e E){i + q e A)} and 

nk-, k ^ + q} r\ k 



i.q.v 



Uk + v, k = ieA (D.5.2) 
rik — V, k = i + q £ A 
Equation (jP.S.ip is very similar to Equation ()D.4.ip . and by repeating the procedure 
used above with Equation (lD.5.ip in place of Equation (]D.4.ip . it can easily be seen 
that a and p must satisfy Equation ()D.4.6P and so that ()D.4.9P is a necessary condition 
for the existence of a product measure steady state of the process on Z'^ as well as on 
a finite torus. 

D.6 Conclusion 

We have shown that there is a straightforward necessary and sufficient condition, Equa- 
tion (|D.1.10p . for a generalized Zero Range Process to have a product measure steady 
state. For Mass Transport Processes, we have found a condition, Equation ()D.2.9p . for 
the existence of a product measure steady state, which is considerably more opaque 
than in the GZRP; it is not clear that this is equivalent to the sufficient condition 
expressed in Equation ()D.2.12p , the counterpart of the condition we have obtained for 
GZRPs. We have, however, presented some reasons to believe that it is. 
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